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SUMMARY 

Ensemble, time, and space averages as applied to turbulent quantities are discussed, and pertinent 
properties of the averages are obtained. Those properties, together with Reynolds decomposition, are 
used to derive the averaged equations of motion and the one- and two-point moment or correlation equa- 
tions. The terms in the various equations are interpreted. The closure problem of the averaged equations 
is discussed, and possible closure schemes are considered. Those schemes usually require an input of sup- 
plemental information unless the averaged equations are closed by calculating their terms by a numerical 
solution of the original unaveraged equations. The law of the wall for velocities and temperatures, the 
velocity- and temperature-defect laws, and the logarithmic laws for velocities and temperatures are derived. 
Various notions of randomness and their relation to turbulence are considered in the light of ergodic theory. 


INTRODUCTION 

Although the unaveraged equations of the last chapter can, in principle, be applied directly to 
turbulence, the historical tendency has been to average out the fluctuations so as to obtain simply vary- 
ing functions. The idea is that simply varying nonrandom values are easier to deal with than the hap- 
hazard motion characteristic of turbulent flow. Physically relevant equations are obtained in this way, 
but the price to be paid as far as obtaining solutions is concerned (the closure problem) is considerable, as 
will be seen later in this chapter. First we will consider various averages of the turbulent quantities. 


4.1 AVERAGE VALUES AND THEIR PROPERTIES 

For the most general turbulent flows an ensemble average over a large number of macroscopically 
(but not microscopically) identical flows is appropriate. In all of those flows, the macroscopic determin- 
ing parameters (e.g., mean velocity, scales, etc.), but not fluctuating quantities, are the same. The 
ensemble average of a quantity, say the velocity, at a point and time t is the arithmetic average of 
that quantity for all the flows. For instance, for a velocity component Uj at point and time t, 

N 

u i( x k> tJn = lim E U i( n ’ x k> t)/ N > ^ 4_1 

n-i ~ n =l 


where n indicates the nth flow, N is the total number of flows over which the average is taken, and the 
overbar with the subscript n designates the ensemble average (the average over n). Needless to say, 
this type of average, although often used in theoretical work, would be hard to implement experimentally, 
because a large number of macroscopically identical flows would not likely be available. 



Fortunately, in most cases statistical uniformity or stationarity with respect to one or more coordi- 
nates and/or with respect to time obtains. Then the average is taken with respect to the one or more 
coordinates and/or with respect to time. For instance, if the turbulence is statistically stationary with 
respect to time (if the ensemble average does not vary with time t), the time average of Uj at a point 
x k is 


u i( x k)t = lim Jo u i( x k> t)dt/T, 

T — > w 


(4-2) 


where the time average is designated by the overbar and the subscript t. 


If the turbulence is statistically stationary with respect to one or more variable coordinates Xj (if 
the ensemble average does not vary with the one or more coordinates Xj), then the space average with 
respect to Xj at a time t and at fixed coordinate(s) x k , where j * k, is 


Uj(x k , ty x . = lim J Xj x Ui(x k , t, Xj )d Xj / (2Xj) , 


(4-3) 


x i 


where the space average is designated by the overbar and the subscript X j, and the integration is over the 
one or more coordinates Xj for which statistical uniformity obtains. (Note that if j can have the values 
1, 2, and 3, x k will be absent from equation (4-3), since k ¥> j.) 

A space average also makes sense for periodic boundary conditions, even if the turbulence is not 
statistically uniform over one or more coordinates. In that case as in the preceding case, the space- 
averaged quantities (averaged over a period) do not vary with position, since it makes no difference where 

the period starts (3uj(x k , t) x / dxj = 0, k * j). Then 

Ui(x k , t)“ = J Xj x Ui (x k , t, Xj)dxj / (2Xj) , (4-4) 


where 2Xj is one spatial period. 

Finally, it might be useful to take the average over both time and one or more coordinates, when 
the turbulence is statically stationary with respect to both. Then the average with respect to t and Xj is 


u i( x k)t, Xj = 


lim J-x Jo T u i< x k' *• X i> dt dX i A 2X i T) • 

Xj,T-> » 1 


(4-5) 


Of course, if the boundary conditions in equation (4-5) are periodic in some of the directions Xj, the aver- 
age in those directions need be taken only over one period, as in equation (4-4). 


The averaging processes considered in this section have been illustrated by taking average values of 
a velocity component Uj. The same processes can, of course, be applied to the mechanical pressure o 
(eq. (3-14)) and to functions of velocities and of pressures such as UjUj, UjUjU k , UjUjU k , and <ru k , where 
the unprimed, primed, and double-primed quantities refer to values at different spatial points. The Uj in 
equations (4-1) to (4-5) need only be replaced by <7, UjUj, etc. It might be mentioned that the simplest 
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average which is descriptive of turbulence is not that of u-, which describes the overall flow rather than 
the turbulence, but that of the second-order tensor UjUj. 


4.1.1 Ergodic Theory and the Randomness of Turbulence 

Turbulence is generally taken to be ergodic, in which case the ensemble, time, and space averages 
of a turbulent quantity (say UjUj) should all have the same value, assuming that those averages exist. 

An ergodic system embodies the weakest notion of randomness in a hierarchy of systems (refs. 1 and 2). 
The so-called mixing systems (those whose variables become uncorrelated as their temporal separation 
At — *• °°) have a stronger notion of randomness than do those that are only ergodic, and systems that 
exhibit sensitive dependence on initial conditions, or chaoticity, have a stronger notion of randomness 
than do those that are only ergodic or only ergodic and mixing. Mixing implies ergodicity, and chaoticity 
implies both ergodicity and mixing, but the converse is not true. At the top of the hierarchy are the 
most random systems; those that, though deterministic, may appear in a certain sense to behave as ran- 
domly as the numbers produced by a roulette wheel (ref. 1) 

So in order for a flow to be identified as being random (or apparently random) in some sense, and 
thus as turbulent, it must be at least ergodic. As mentioned in chapter I, turbulent systems appear to be 
at least as random as chaotic systems. In chapter VI it is argued that they are likely more random. 

Because of the ergodicity of turbulence a distincti on among the various kinds of averages will not 
usually be made, and they will be written simply as Uj, UjUj, U;<7' , etc., the subscripts n, t and x- in 
equations (4-1) to (4-5) being omitted. It will always be assumed, of course, that the averages taken are 
of an appropriate type, in line with equations (4-1) to (4-5). 


4.1.2 Remarks 

The averaging considered here, which is known as Reynolds averaging, is the type which will be 
used in this book. There are, of course, other types of averaging, which have specific uses. For instance 
conditional averaging, in which averages are taken under some specified condition, such as the condition 
that only velocities greater than some value be used in the average, is sometimes useful. However none of 
the methods of averaging circumvents the closure problem (considered in section 4.3) which occurs when 
the nonlinear continuum equations are averaged. 

4.1.3 Properties of Averaged Values 

Finally, we will consider some properties of averaged values which will be useful in obtaining the 
averaged equations of fluid motion. We note first that the derivative of an average equals the average of 

the derivative. Thus, for example, dujUj/8x k = 3ujUj/dx k . This equation can be obtained from equa- 
tion (4-2) (with Uj replaced by UjUj) if the average is taken over time. Thus 


3 



(4-6) 


d“i u j 

dx t 


dx 


kj 


, Um Jo ^ U i U j) dt / T 


T — > w 


i im lc 


|*T 

d(u iUj ) 

Jo 

dx k 


dt/T 




dx b 


The same result is obtained if the average is other than that over time. Also, it is easy to show that the 
sum of averages equals the average of the sum. Thus, again using equation (4-2), 


UjUj' + u(u^ = lim (ujUj)dt /T + lim J 1 (u/u^dt/T 
T— > PO T-+ OO 

= lim Jq (uju/ + Uj Uj )dt /T = Ujuj + u-uj . 


T-> <*> 


(4-7) 


Moreover, taking an average of an average (designated by a double bar) does not change its value, since, 
for example, 


UjUj 


lim (ujUj)dt / T = UjUj lim dt/T = UjUj . 


T — > w 


T— > 


(4-8) 


Finally, the average of the product of an averaged and an unaveraged quantity is the product of the first 
(averaged) quantity by the second quantity averaged. For example, 


u i( u j u k) = lim J J Ui(u k Uj)dt / T = uj lim F (u k Uj)dt / T 

T— > « T— ► » 


(4-9) 


u i u j u k 


As an example of a relation for averages which does not hold, we note that 


since 



All of the preceding relations which were shown to hold for time averages hold, of course, for the 
other kinds of averages. Thus, for instance, if we consider ensemble averages and use equation (4-1), we 
find that 
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(4-8a) 


N 

U|Uj = lim £ ujuj/N = ujm, 

N- « n=l 

which is the same result as that obtained in equation (4-8) for time averages. 


4.2 EQUATIONS IN TERMS OF MEAN AND FLUCTUATING COMPONENTS 


Equations for mass, momentum, and energy conservation, as well as for the pressure, which are 
general enough for all of the work in this book are given by equations (3-4), (3-18a), (3-34), and (3-21a) 
respectively. We first rewrite those equations in slightly different forms and notations: 


JP| 

II 

o 

(4-10) 

*IJ< 

II 

1 

d( 3 i S k) 1 ~ <r e ) dV 

- - 3- + ", ' - HI ~ T.) gi , 

«x k p dxj dx k ox k 

(4-11) 


8T 5(Tu k ) 

= - + a , 

dt 3x k <3x k dx k 

(4-12) 

and 

1 - * e ) _ _ ^( Q i Q k) _ ^ cT 

p dxj dx ( 9xj 5x k 1 dxj 

(4-13) 

These equations are the same as equations (3-4), (3-18a), (3-34), and (3-21a), except that ~’s have been 
placed over the instantaneous velocities Uj, mechanical pressures a, and temperatures T, and the non- 
linear terms have been written in the so-called conservative form by using the continuity equation (4-10). 

Following Reynolds (ref. 3), one can break the instantaneous quantities in equations (4-10) to (4-13) 
into mean and fluctuating (or turbulent) components. That process is known as Reynolds decomposition. 
Thus, set 


Ui = Uj + u; 

(4-14) 


a = P + a 

(4-15) 

and 

T = T + r 

(4-16) 
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where the first and second terms on the right sides of the equations are respectively mean and fluctuating 
components, and where 



(4-17) 


U, 




(4-18) 


F = 7, (4-19) 

and 

T — T. (4-20) 

As usual, the overbars designate averaged values. Equation (4-10) becomes, on using equations (4-14), 
(4-17), (4-18), and the properties of averaged values given in equations (4-6) to (4-8), 


6U k = _ n 
dx k dx k 


(4-21) 


which shows that both the mean and fluctuating velocity components satisfy conservation of mass. Equa- 
tions (4-11) to (4-13) become, on using equations (4-14) to (4-20) and (4-6) to (4-8), taking averages, and 
subtracting the averaged equations from the unaveraged ones, 


duj 

at 


9 ( u i u k) 1 d(<r - <r e ) 


3Ky 


fa; 


# 11 : dUi dU: 

+ U - /3gj7- -u k - U k + 

dx k dx k dx k dx k dx k 


dUiUfc 


(4-22) 


and 


dr 

at 


5 (™k) ^ a 2 r 

ax k dx k dx k 


u k ^-u k ^ + !^, 
dx k ax k ax k 


(4-23) 


i a(<r - <r e ) _ _ ^(uj^) dr auj au k &up k 

' ' " " — ~j~ — • 

p dxj dxj dxj dx k dxj 3x k dxj dxj 5x k 


(4-24) 


Equations (4-22) to (4-24) are useful for studying turbulence processes and for constructing evolution 
equations for correlations (e.g., for UjUj). The first five terms of equation (4-22), the first three of equa- 
tion (4-23), and the first three of equation (4-24) are similar to the terms in equations (4-11) to (4-13) 
respectively, although their meanings are exactly the same only if U- = T = P = 0 (see eqs. (4-14) 
to (4-16)). 

As was the case for equation (3-19) (or (4-11)), the first three terms on the right side of equation 
(4-22), which give contributions to the rate of change of the velocity fluctuation dUj/at, can be inter- 
preted as an inertia-force (or turbulence self-interaction) term, a pressure-force term, and a viscous-force 
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term. The remaining terms are, respectively, a buoyancy-force term, a turbulence-production term, a 
mean-flow convection term, and a mean turbulent stress term, which may appear when the turbulence is 

statistically inhomogeneous (when mean turbulence quantities such as UjU k are functions of position). 
(The reasons for referring to the production and convection terms as such will perhaps become clearer 
when the equivalent terms in the averaged equations are discussed.) It will be seen in the next chapter 
that when the mean-velocity gradient is not zero, the term — Ujjduj/dx^ generates a small-scale struc- 
ture in the turbulence by vortex stretching or by a breakup of eddies into smaller ones. The nonlinear 
self-interaction term — d(UjU k )/dx k also produces a small-scale structure, and in addition produces ran- 
domization of the flow. This effect will also be considered in the next chapter. 

The terms on the right side of equation (4-23), which give contributions to the rate of change of 
the temperature fluctuation dr/d t, are respectively a nonlinear temperature-velocity interaction term, a 
molecular diffusion term for temperature fluctuations, a production term for temperature fluctuations, a 
mean-flow convection term for temperature fluctuations, and a mean turbulent heat-transfer term which 
may be present when the turbulence is inhomogeneous. 

The Poisson equation for the mechanical pressure fluctuations has four source terms — a nonlinear 
term, a buoyancy term, a mean-velocity-gradient term, and a mean-turbulent-stress term which appears 
when the turbulence is inhomogeneous. 

4.3 AVERAGED EQUATIONS 

Although the averaged equations for fluid motion do not form a closed (complete) set, they are 
very useful for studying the physical processes in turbulence. Moreover, approximate solutions can be 
obtained by introducing various closure schemes. Alternatively (or preferably), the terms in the averaged 
equations can often be calculated from a numerical solution of the unaveraged equations. 


4.3.1 Equations for Mean Flow and Mean Temperature 

First consider the equations obtained by averaging each term in equations (4-11) to (4-13) after 
applying Reynolds decomposition (eqs. (4-14) to (4-20)) and using (4-21). This gives 


auj 

at 


/>u k 


auj a(p - » e ) 


ax b 


ax ; 


a 


f au i , 

au k 


+ — 

pv 


dX: 

1 V 

- /> U i U k 

ax k 

dx k 


-^(T-T.)g,, (4-25) 


aT TT aT a 

pc— = — />cU k 


at 


a*, ax b 


t \ 

aT , — 

— pea + /»cu k r 

ax k 


(4-26) 


and 


d*(P ~ ».) 

d)(j dx g 


eJ a 

“ ~~ 


u b 


au ( 


dx 


kj 


a% u k a ar 

" r- - 

OX | ax k ox • 


(4-27) 
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where the properties of averages given by equations (4-6) to (4-9) were used. One should note that time 
averages would not be appropriate in these equations unless averaged quantities change much more 
slowly with respect to time (preferably infinitely more slowly) than the unaveraged quantities. Other- 
wise, averages with respect to space variables in directions for which the turbulence is statistically sta- 
tionary would be preferable. Equation (4-26) is written for a liquid; for a perfect gas c should be 
replaced by c p (see equation (3-35)). Equations (4-25) to (4-27) look like (4-11) to (4-13) with instan- 
taneous values replaced by average values, but with the important difference that an extra term involving 
UjU k or u k r now appears in each of the equations. These terms arise from the nonlinear velocity and 
velocity-temperature terms in equations (4-11) to (4-13) and are a manifestation of the closure problem in 
turbulence. If those terms were absent, the set of equations (4-25) to (4-27) would contain as many 
unknowns as equations and so could be solved. In that case turbulent flows would be no more difficult to 
calculate than laminar ones. Note that barred quantities in equations (4-25) to (4-27) which contain 
lower-cas§ letters are turbulent quantities. The term — puju k was discovered by O. Reynolds (ref. 3), 
aftd is often called the Reynolds term or the Reynolds stress. 

4. 3. 1.1 Interpretation of the terms — />UjU k and />cu k r. — The forms of equations (4-25) and (4-26) 
suggest that the quantities — /?UjU k and pcu k r respectively augment the molecular (or viscous) stress 
tensor pi/(dU- l /dx k + dU^/cbcj) and the molecular heat transfer vector — pcacfT/dx^. Since they involve 
fluctuating or turbulent components, we interpret — /?UjU k as a turbulent (or Reynolds) stress tensor and 
pcu k r as a turbulent heat transfer vector. (The quantity UjU k is a second-order tensor, since it is the 
average value of the product of two vectors (the average value of a tensor is a tensor), and u k r is a vec- 
tor, since it is the average value of the product of a vector and a scalar (the average value of a vector is a 
vector)). Thus in equations (4-25) and (4-26) we write 


and 



fdUj 

dtV 

pv 

— + 


ax k 

axjJ 


/> u i u k 


(4-25a) 




Q k = -pea + pcu k r, 

dx k 


(4-26a) 


where r ki is the total stress tensor (the sum of the molecular and turbulent stress tensors), and Q k is 
the total heat transfer vector (the sum of the molecular and turbulent heat transfer vectors). 

We note that the expression for the turbulent stress corresponds exactly with that for the mole- 
cular stress obtained in the kinetic theory of gases (see, e.g., ref. 4). It is only necessary to replace the 

macroscopic turbulent velocity fluctuations in — /?UjU k by random molecular velocities. A similar corre- 
spondence exists in the expressions for the turbulent and molecular heat transfer, where the macroscopic 
velocity and temperature fluctuations in /£>cu k r are respectively replaced by a random molecular velocity 
and molecular kinetic energy (molecular temperature) (ref. 4). 

Consider now the term — />UjU k for i * k. Then, for instance, — iffu 1 u 2 will, in the presence of a 
mean-velocity gradient dUj/dxj, act like a turbulent shear stress on an Xj - x 3 plane. Similarly, /»cu 2 r 
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will, in the presence of a mean-temperature gradient 3T/dx 2 , act like a turbulent heat transfer in the 
— x 2 -direction. 

To see how these effects come about, consider figure 4-1, where the curve represents either Uj or T 
plotted against x 2 . Then if u 2 is positive at a particular location, the fluid will instantaneously be moving 
into regions of higher U 1 (if dU 1 /dx 2 is positive as shown), and so u x will tend to be negative (the local 
velocity will tend to be less than UJ. Similarly, if u 2 is negative, the fluid will, instantaneously, be moving 
into regions of lower U 1; and so Uj will tend to be positive. In both cases the product u^ will tend 

to be negative, and so the average value u 1 u 2 will be nonzero and negative. That is, there will be a 
negative correlation between u^ and u 2 . The quantity — pu^u 2 is the turbulent (Reynolds) shear stress 
acting on an Xj — x 3 plane and augments the molecular or viscous shear stress pi/dU|,/dx 2 on that plane. 
Of course if Uj is uniform, UjU 2 will be zero; there will be no correlation between u x and u 2 . 

A similar interpretation applies to pcn^r if the curve in figure 4-1 represents the mean tempera- 
ture T plotted against x 2 . That is, when u 2 is positive at a particular Xj, the fluid will, instantane- 
ously, be moving into regions of higher T (if dT/dx 2 is positive as shown), and so r will tend to be 
negative (the local temperature will tend to be less than T). Similarly, when u 2 is negative t will 

tend to be positive. In both cases the product u 2 r will tend to be negative. The average value u 2 r 
will then be nonzero and negative, so that like the molecular conduction term, pc\i 2 r will produce heat 
transfer in the -x 2 -direction. If T is uniform there will, of course, be no correlation between u 2 and 

r, and u 2 r will be zero. 


Finally, we should say something about — /JUjU^ when i = k. In that case — puju^jj (no sum on 
i) will act like a normal turbulent stress and will augment the normal molecular stress pvdU-Jdx^y For 
2 

instance -pu 1 will be the normal stress on an x 2 — x 3 plane. 

The existence of the turbulent stress tensor -pUjUj and heat transfer vector pcruj is an impor- 
tant and physically significant deduction from the continuum equations for fluids (the Navier-Stokes, 
energy, and continuity equations). The deduction was obtained from those equations with no approxima- 
tions by using only Reynolds decomposition and the rules for averaging. Unfortunately the procedure 
does not provide a way of calculating the values of the turbulent stresses and heat transfer, and so we are 
left with a closure problem (more unknowns than equations). Next we will consider some simple closure 
schemes. 


4.3.2 Simple Closures of the Equations for Mean Flow and Temperature 

In order to close the system of equations (4-25) to (4-27), that is, to write it in a form in whi ch the 
number of unknowns equals the number of equations, one must write the quantities UjUj and u^ r in 
terms of the mean velocity Uj, the mean temperature T, and Xj. To do that we necessarily introduce 
additional information into the equations, so that the theory is not deductive. However, because of the 
practical importance of obtaining solutions, a great many closure proposals have been made. 
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Probably the simplest way of closing the system of averaged equations would be to assume that ujuT 
and u k r are respectively proportional to UjUj and U k T. One might suppose that this is a reasonable 
assumption because UjUj and U k T have the same dimensions as UjUj and u k r and are respectively 
a second-order tensor and a vector (as are UjUj and u)~r). But consider what happens to the turbulent 

stress terms in equations (4-25) and (4-27) for the simple case in which Uj = tf.jUj, and U 1 is independent 
of Xj; that is, the flow is fully developed. In that case the suggested closure assumption gives, in equa- 
tions (4-25) and (4-27), d(^UjU k )/ dx k = d^UjUj)/ dx 1 = 0, and d 2 ^^/ dx ( dx k = ^(UjUj)/ dx^ = 0, 
where p is constant. Thus the assumed form for UjUj would not appear in equations (4-25) and (4-27), 
and there would be no effect of turbulence on fully developed flow. Since that result is contrary to exper- 
ience, the proposed closure assumption cannot be reasonable. Evidently the conditions of correct dimen- 
sionality and correct tensor properties, while necessary in a closure assumption, are not sufficient. 


4.3.2. 1 Eddy diffusivities .- — Inadequacies such as that noted in the expression for UjU: which was 
just considered can be avoided by introducing the so-called eddy diffusivities. This can be done formally 
by replacing molecular quantities in equations (3-13), (3-14), and (3-29) by turbulent quantities. Thus, 

in those equations, <Tjj / p —* — UjUj, p/p — > e, qj/(pc) — > ruj, and k/(pc) — > where the arrows are 

read “becomes.” In addition, according to our notation for turbulent flow, the mean velocity is written 

as Uj. Then, for incompressible flow we get, for — UjUj and ru7, 


and 


UjUj - * 


(dUj aUj 

dxj dxj 


- u k u k ^ij 


— _ dT 

ru j - — £ h -jr— > 

dxj 


(4-28) 


(4-29) 


where e is variously called the turbulent viscosity, eddy viscosity, or eddy diffusivity for momentum 
transfer, and is the turbulent conductivity, eddy conductivity, or eddy diffusivity for heat transfer. 

Unlike the molecular viscosity and conductivity in equations (3-13), (3-14), and (3-29), the eddy 
diffusivities are functions of the character and intensity of the turbulence. We notice that equations (4-28) 
and (4-29) do not suffer from the problems associated with the expressions given near the beginning of sec- 
tion 4.3.2; in general they give nonzero values for both the turbulent stress and heat-transfer terms in 
equations (4-25) to (4-27) for a fully developed flow in a channel or pipe. Moreover they give zero values 
for UjU 2 and ru 2 when the mean velocity and temperature gradients are zero, as they should at 
the center of symmetrical flows (see section 4.3. 1.1). For unsymmetrical flows equations (4-28) and (4-29) 
may break down, but we shall not be concerned here with those flows. For statistically homogeneous 
flows without mean velocity gradients equation (4-28) gives UjUj = u k u k 6y/3 , which is true only for iso- 
tropic turbulence. However for homogeneous turbulence we will use other methods. The eddy diffusivi- 
ties in equations (4-28) and (4-29) are scalars. Other expressions, in which they are tensors, have also 
been proposed (ref. 5). 


Equations (4-28) and (4-29) do not provide closures for the equations (4-25) to (4-27), since we still 
do not know e and as functions of position and/or mean velocities. However they provide a framework 
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in which UjUj and raj are respectively a second-order tensor and a vector, as they should be, and in 
which reasonable expressions for UjUj and ruj might be obtained. In general, expressions for those 
quantities must be tailored to the particular problem being considered. 

4.3. 2. 2 Mixing length . — Before considering specific expressions for the eddy diffusivities, we will intro- 
duce Prandtl’s mixing-length hypothesis (refs. 6 and 5) which gives a rough estimate of e^/e and which 
gives an approximate picture of how momentum and heat might be transferred in a turbulent shear flow. 

When turbulence exists in a flow, eddies or portions of fluid move about in an apparently random 
fashion. If a mean velocity gradient and/or temperature gradient exist in a direction transverse to the 
main flow, some of the eddies will move transversely into layers of different mean velocity and/or temper- 
ature (see fig. 4-1 and sec. 4. 3. 1.1). Consider the turbulent fluctuations u 1? and r at a point x ; (see 
eqs. (4-14) and (4-16). The mean velocity and temperature at that point are respectively Uj and T, 
and the mean gradients are in the x 2 -direction, as in figure 4-1. 

According to mixing-length theory the fluctuations u x and t at Xj are produced by an eddy or 
turbulent particle which originates from an instability at another point. We define the virtual origin of 
the eddy that produced the fluctuations Uj and r at Xj as the point Xj 0 where the eddy would have 
been born with the mean velocity Uj 0 and temperature T 0 if the Xj-moraentum and temperature of 
the eddy were conserved as it travels from Xj 0 to Xj (see fig. 4-1). Note that actual conservation need 
not occur, since the virtual origin of the eddy can differ from its actual origin. 

From equations (4-14) and (4-16), which define Reynolds decomposition, the fluctuations u x and 
t at Xj are given by 


ui = ui - Uj 

and 

r =T - T, 


(4-30) 

(4-31) 


where the ~’s designate total instantaneous quantities at Xj, or mean values plus fluctuations. But mean 
values plus fluctuations in equations (4-30) and (4-31), designated by Uj and T respectively, are just 
the total Xj-velocity and temperature of the eddy considered in figure 4-1 and in the definition of virtual 
origin in the last paragraph. According to that definition, if the eddy were born at the virtual origin Xj 0 
with the mean Xj-velocity Uj 0 and temperature T Q , then the Xj-momentum and temperature of that 
eddy would be effectively conserved as it moves from Xj 0 to Xj. Thus Uj and T can be replaced by 
U 10 and T 0 respectively. Equations (4-30) and (4-31)’ then become 

(4-32) 
(4-33) 


and 


*1 — u 1)0 


U, 


T 0 - T. 


Substituting equations (4-32) and (4-33) into -u^ and ru 2 , we get 

~ u l u 2 = (Ui~ U 10 )u 2 


(4-34) 
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and 


m 2 — (T Tq) u 2 • 


(4-35) 


Note that Uj — Uj 0 = Uj and T — T 0 = t are not removed from the bar since they vary with the 
averaging variable(s’) (xj, x 3 , and/or t). If we expand Uj and T in Taylor series about Uj 0 and T 0 
respectively and retain only the first two terms in each, we get 


and 


dUj dU, 

~ u l u 2 = *2 — »2 = « 2 U 2 — 
CIX2 QX2 


(4-36) 


ru 2 



M 


2 U 2 


dT 

dx 2 


(4-37) 


where ( 2 = x 2 - x 2 0 and is designated the mixing length. Let i = 1, j =2, and the mean flow be in 
the Xj-direction in equations (4-28) and (4-29). Comparison of those equations with equations (4-36) 
and (4-37) then gives 


e = = «2 U 2 • 


(4-38) 


The relation between e and given by equation (4-38), besides following from mixing length theory, 
gives good agreement with experiment except at low Prandtl numbers cju/k. At low Prandtl numbers, 
as for liquid metals, the thermal conductivity k is generally so high that heat is conducted to or from an 
eddy as it moves transversely. The temperatures in equation (4-35) will then have to be replaced by val- 
ues which are closer together, and | ru 2 | will decrease (refs. 7 to 9). In this chapter we will not be con- 
cerned with low Prandtl number fluids and will use equation (4-38). 

Equations (4-36) and (4-37) have sometimes been criticized (often mercilessly) because they assume 
that f 2 is small enough that Uj and T vary linearly over that distance, whereas in reality it may not 
be. That problem could be overcome by retaining more terms in the Taylor-series expansions for 
Uj — U 10 and T — T Q . For instance in the expansion for Uj — Uj 0 , we could retain an additional 
term (l/2)^6 2 U 1 /dx 2 . However that would complicate the analysis and may not be worth the effort, 
since one can probably absorb any second-order effects in the expressions which are assumed for the eddy 
diffusivity. 


4.3.2. 3 The nonuniformity of turbulent mixing . — A fundamental question about the nature of 
turbulence concerns how turbulent mixing takes place. Here we consider the instantaneous turbulent 
mixing which occurs in the presence of mean velocity and/or temperature gradients. We consider it at 
this point because of its relevance to mixing-length theory. It turns out that nonuniform mixing is a 
consequence of several known facts about turbulence. 

Let us see what the presence of the turbulent stress and heat transfer terms in the equations for 
mean flow and heat transfer (eqs. (4-25) and (4-26)) implies about the instantaneous turbulent mixing. 
Instantaneous mixing refers here to the mixing one would see in a snapshot taken at a particular time. 
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Note first what would happen if the spatial pattern of instantaneous turbulent mixing were uni- 
form, or nearly so. If that were the case, a portion of fluid as it moves transversely in mean velocity and 
temperature gradients (mean velocity in the Xj-direction) would have a uniform tendency (because of 
uniform mixing) to assume the mean Xj-momentum and temperature of the surrounding fluid at each 
point along its path. That tendency would be more pronounced at higher turbulence intensity or Reynolds 
number, because small-scale motions become excited with increasing turbulence intensity (ref. 10), and so 
the turbulent mixing (average or instantaneous) increases. (Note that turbulent mixing takes place most 
efficiently by small-scale motions, since those provide the most intimate contact of the fluid entering a 
region with that already there.) 

Thus if the instantaneous turbulent mixing were spatially uniform, the tendency of a portion of 
fluid to assume the mean Xj-momentum and temperature of the surrounding fluid at each point as it 
moves transversely would increase with increasing turbulence intensity. That would, however, cause the 
fluctuations from the mean, u 1 and r in the turbulent stress - pu jUj and heat transfer pc ruj , to 
decrease in magnitude with increasing turbulence intensity or Reynolds number. The stress component 

pu x would then decrease. But that trend is unphysical and does not occur. In fact, as might be expected, 
the opposite trend occurs; as turbulence intensity (uju^/3) 1 / 2 or Reynolds number increases, u^ increases 
(ref. 11). 

The instantaneous turbulent mixing therefore cannot be spatially uniform, or nearly so, as assumed 
in obtaining the above unphysical trend. There must be regions of relative quiescence if x 1 -momentum 
and heat are to be transferred turbulently at high turbulence intensities. But in that case there must also 
be regions where the instantaneous mixing is relatively intense and localized, since that is the only way 
the average mixing could be high for high turbulence intensities when regions of quiescence are present. 

So the only sensible assumption about the instantaneous mixing is that it is small except in localized 
regions, where it is intense. Then the above unphysical trends do not occur, since the tendency of a 
portion of fluid, as it travels transversely, to assume the mean x 1 -momentum and temperature of the 
surrounding fluid is sudden and is confined to localized regions. Note that even if the turbulence is 
statistically homogeneous, the instantaneous turbulent mixing tends to be highly inhomogeneous. 

The fact that instantaneous turbulent mixing takes place mainly in localized regions means that 
the turbulence must be spatially intermittent in the small scales, since mixing, as mentioned before, takes 
place mainly by small-scale motions. Intermittency in the small scales has been found experimentally 
(ref. 12). 

The localness or suddenness of the turbulent transfer considered in this section also seems to be in 
agreement with the concept of bursting coherent structures in shear flow near a wall. Much work has 
recently been done on that phenomenon (ref. 13). 

Our result that instantaneous turbulent mixing is sudden and localized is congruous with mixing- 
length theory, which requires a certain suddenness in the turbulent mixing for turbulent transfer to take 
place. The mixing length can be thought of as the effective distance an eddy moves before mixing with 
the surrounding fluid. If the mixing took place continuously the mixing length would be zero, and the 
turbulent shear stress and heat transfer would be zero (eqs. (4-36) and (4-37)). 

Thus, although fluid turbulence occurs in a continuum, changes in the momentum and temperature 
of a moving portion of fluid tend to be sudden and localized. In that respect turbulent systems are not 
unlike the systems considered in the kinetic theory of gases, where encounters between particles are 
sudden and localized. It is of interest that it was apparently kinetic theory that originally inspired 
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turbulent mixing-length theory; the mixing length was supposed to be something like the mean free path 
of kinetic theory (ref. 6). 

4.3. 2. 4 Some conditions satisfied by the turbulent shear stress and heat transfer near a wall. — It 
may be helpful in obtaining models for the turbulent shear stress and heat transfer (or for the eddy 
diffusivities for momentum and heat transfer) to determine conditions which must be satisfied by those 
quantities in the region near a wall. Moreover the results obtained in this section, so far as they go, are 
exact, depending only on the continuity equation and the boundary conditions for the velocity and tem- 
perature fluctuations. The results for the turbulent shear stress have been given previously, for example 
in reference 14. 

Consider an incompressible nonrarefied flow in the Xj-direction parallel to a wall. The normal to 
the wall is in the X 2 -direction, with x 2 = 0 at the wall. The continuity equation is 

duj duj^ du 2 du 3 

dxj dx 1 dx 2 dx 3 

Since the flow is nonrarefied, the nonslip condition holds at the wall, so that 

( u l)x 2 =0 = ( u 2)x 2 =0 = ( u 3)x 2 =0 = 0 (4-40) 


(4-39) 


at every point on the wall. Thus 


duj 

dx-, 


/x 2 =0 




0 


= 0 . 


Equation (4-39) then gives, at the wall, 




=0 


= 0 . 


(4-41) 


(4-42) 


Consider now the turbulent shear stress -/m 1 u 2 . Then 


^l u 2 

dx 2 


dui u 2 

dx 2 



9uj 



(4-43) 


3\u 2 


dx n 


^ u 2 _ ^1 ^2 
— u 1 — 2 

dxl ^2 ^2 

L 


d 2 ^ 



(4-44) 
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and 


d\u 2 

dx 2 


= -u, 


d 3 u. 


dx„ 


diii d 2 u 2 
^2 fo 2 


du 2 d 2 u 1 d 3 Uj 

3 u,, 


^2 ax 2 


dx 


(4-45) 


where the properties of averages given by equations (4-6) and (4-7) were used. By using equations (4-40), 
(4-42), and (4-43) to (4-45) we get, at the wall, 


-( u l u 2)x 3 =o = 0, 


du^ 

dx 0 


= 0, 


' x ,=0 


(4-46) 

(4-47) 


and 


d 2 u 1 u 2 


dx„ 


0, 


/x 2 =0 


fd^uiujl 


= -3 


duj dVjj 


dx„ 


/x-=o 


dx, dx' 


^ Jx 0 — 0 


(4-48) 


(4-49) 


Since the zero’th, first, and second derivatives are zero at the wall, — Uiu 2 , starting at zero, increases 
very slowly with x 2 near the wall. The third derivative (see eq. (4-49)) gives no information, since at 

least at this stage, we do not know how to evaluate - ( duj / dx 2 d 2 ^ / dx 2 ) . Thus 


0 < 


d 3 UjU 2 


dx 


< 00 . 


A similar result holds for higher-order derivatives. 
As an example, determine whether 


- u iu 2 


OC X 


18/5 

2 


(4-50) 


(4-51) 
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is in agreement with equations (4-46) to (4-50). Equation (4-51) and its derivatives give, at the wall, 

~ ( u l u 2)x 2 =0 = — ( ^ u i u 2 / ^ jc 2)x 2 =0 ~ — ( 5^ u i u 2 / ® x 2^* a =0 = 0 ancl (^UjUj/ dx^-Q = <*>. Since these 
evaluations are in agreement with equations (4-46) to (4-50), equation (4-51) satisfies continuity and the 
nonslip boundary conditions for Uj at the wall. In section 4.3. 2.9 we will see that it also gives results in 
agreement with experiment. 

Next consider the turbulent heat transfer pcru^. Although the procedure is similar to that for the 
turbulent shear stress, the results can be different because the temperature fluctuation r at the wall may 
not be zero if the thermal conductivity of the wall material is not sufficiently high. Thus consider two cases, 

r x 2 =0 * 0 (4-52) 

for low thermal conductivity of the wall material, and 

r o = 0 (4-53) 


for high thermal conductivity of the wall material. Then, proceeding as for u ^u 2 , we have 


dru, drun 


dx 0 


dxn 


= r . 


du 

dx 


2 . dr 

+ u 2 

dXn 


(4-54) 


d 2 




d 2 u 


dx! 


= r. 


2 + 2 dT ^ 


+ u 

dx\ ^ 2^2 


and 


d 3 


TU„ 


ax 


= + 3^1 fl + u.£l 

K 8X2 *2 8X2 *2 


(4-55) 


(4-56) 


Using equations (4-40), (4-42), (4-52), and (4-54) to (4-56), we get, for a wall with low thermal 
conductivity, 


(™2)x,=0 = 0 , 


(dru 2 
dx 2 

v 1 J 


= o, 


r x 2 =0 


(4-57) 


(4-58) 
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(4-59) 


2 V° 


dr ® u 2 

^2 ax 2 

^ 1 Jx 2 =0 


Again, since we cannot evaluate the right side of equation (4-59), 


2 L _n 


A similar result holds for the third and higher derivatives. 

On the other hand, equations (4-40), (4-42), (4-53), and (4-54) to (4-56) give, for a wall with high 
thermal conductivity, 

~ (™ 2 ) x 2 =0 = 0 > ( 4 ‘ 61 


^2 Ac I 
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and 


0 < 


d 3 ? 


& 


< 00 


bc 2 =° 


(4-65) 


A result similar to equation (4-65) holds for the fourth and higher derivatives. 

So, like the turbulent shear stress, the turbulent heat transfer is zero at a wall and increases very 
slowly with x 2 near a wall. However |ru 2 | may increase slightly more rapidly than ] u 1 u 2 | if the 
wall thermal conductivity is low, since in that case (crru 2 / 3x 2 ) X2 _ 0 may be nonzero. But since an 
effect of wall thermal conductivity on turbulent heat transfer does not seem to have been measured, and 
since equality of eddy diffusivities for heat and momentum transfer gives good results, this possible dif- 
ference may not be significant in most situations. Of course, new work may show otherwise. 

4.3.2.5 Specialization of equations (4-25) and (4-26) for fully developed parallel mean flow and fully 
developed heat transfer without buoyancy. — For simplicity consider a turbulent flow between two infinite 
parallel walls, which may or may not be moving. The flow is fully developed in space and time and is in 
the Xj-direction. The direction normal to the walls is x 2 . Fully developed flow is here taken to mean 
that all of the dependent variables in equation (4-25) except the pressure are independent of time and Xj. 
The pressure may vary with Xj. If we neglect buoyancy effects by letting gj = 0, equation (4-25) be- 
comes, for i = 1 and 2, 


and 


dP 

dx t 


+ 


dXn 


pv. 


dUj 

dx 0 


PUjUj 


(4-66) 



dx 2 dx 2 


(4-67) 


Note that by virtue of the equation following (3-22), o & drops out of equation (4-25) for gj = 0. Differ- 
entiating equation (4-67) with respect to x x gives 


d 

f \ 

dP 

d 

( 

du 2 





— — P 


dx 2 

dxi 

dx x 

dx 9 
k L ) 


(4-68) 


2 2 

since u 2 = u 2 (x 2 ) does not change with x x for fully developed flow. Integrating equation (4-66) with 
respect to x 2 (using eqs. (4-68) and (4-25a)) then gives 


dP 

dx x 


x 2 


+ C x 


dUj _ 

= pv— — - pu x u 2 = r 21 , 
dx 2 


(4-69) 
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where C x is a constant of integration whose value depends on the boundary conditions, and r 2 1 the 
total averaged shear stress (the sum of the molecular and turbulent parts). According to equation (4-69) 
the total shear stress varies linearly with x 2 . 


For the case where both walls are stationary, dU 1 /dx 2 = u 1 u 2 = r 21 = 0 at the channel center 
(at x 2 = (x 2 ) c ), and 

c, = -(*,), |^ ■ (4-70) 

Using that value for C^, equation (4-69) becomes, at x 2 = 0, 


and 


- ( x 2)c— — ~ ( r 2l)w — r v, 
(TXi 


x 2 _ pu *^1 
( x 2)c r w dx 2 



l u 2 


'21 


(4-71) 


Thus the total shear stress r 21 varies linearly with distance from the wall, being a maximum at the wall 
and zero at the channel center. Since the turbulent shear stress — pu 1 u 2 is much greater than the mole- 
cular shear stress /»i/dU 1 /dx 2 , except in the immediate vicinity of the wall, — UjU 2 also varies linearly 
with x 2 over most of the channel cross-section. At the wall, of course, unlike the total shear stress 

r 2l > ~ u i u 2 i s zero because the velocity fluctuations are zero. 


Next consider the energy equation (eq. (4-26)) for the case to which equation (4-69) applies (a 
steady, turbulent, unidirectional, fully developed mean flow between infinite parallel walls). The mean 
temperature is also steady, and the flow is fully developed thermally. Equation (4-26) can then be 
written using equation (4-26a), as 


/jcUi 


dT 


d 

dT 

\ 

, d 

\ 

dT — 


pea 

pen 

+ — 

pc a — pcu 2 r 

dxj 

dx, 

V. 1 

, 

dx 2 

dx 2 

v £ / 


dQi 

dx. 


dQ 2 

dx 9 


(4-72) 


where Q 1 and Q 2 are respectively the total (time-averaged) heat transfer (molecular plus turbulent) in 
the Xj- and x 2 -directions. Fully developed thermally is here taken to mean that the molecular and tur- 
bulent heat-transfer vectors are independent of Xj. This definition is consistent with equation (4-72), as 
can be seen by differentiation of that equation with respect to Xj. (Note that although the left side of 
equation (4-72) was not assumed independent of Xj, that independence follows from the independence of 
Xj ■ of the molecular heat-transfer vector.) Then equation (4-72) becomes, for fully developed heat transfer, 
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or 


pcU! 


dT 

dx 1 


d 

dx 2 


|pc a 


dT 

dx 2 




- pcn 2 r 



(4-73) 


f U l( x 2) dx 2 + C 2 = <*— “ u 2 r = 

UXj dx 2 



(4-74) 


where dT/dxj is written outside the integral sign, since the flow is fully developed thermally 
((9/^HadT/dXj) = 0 = (d/dx^dT/dxJ). 

4. 3. 2. 6 Law of the wall. — Equations (4-11) and (4-13) can be written in dimensionless form for 
gj = 0 as 


dm __ d(V fi k + ) cb + , ^ G i + 


dt + 


dx, 


dx. 


dx, + dx. + 
k k 


(4-75) 


and 


2~ + 


dV 


dx+dx + 
t I! 


d^iX) 


dx* dx? 

i k 


(4-76) 


where over a quantity indicates a total instantaneous value (mean plus fluctuating), and where 

+ x i f T Jp 

v rr * 


u. 

i 


+ _ u i 


vJp 


- + _ o j p , + tT w 

o - — _, and t = 


T w/P 


(4-77) 


pv 


(The quantity <r e drops out of the set of equations since g; = 0 (see equation following (3-22)).) 

The set of equations (4-75) and (4-76) contains four equations in the four unknowns uj, u' 2 , u.^, 
and a . Thus one could do a thought experiment, starting from given initial conditions. Using appro- 
priate boundary conditions at the walls (e.g., ut = 0 for x 2 = 0, 2(x 2 ) c , a solution of equations (4-75) 
and (4-76) could, in principle, be obtained. For fully developed statistically steady flow one would have 
to consider a region far from the channel entrance and times large enough for the solution to become inde- 
pendent of initial conditions. The parameter in this system is (x 2 ) c , which appears as 1/2 the maximum 

value of x 2 (th e wall s are 2(x 2 ) c units apart). The quantity (x 2 ) c = (x 2 ) c \j r w / pj v is a kind of Reynolds 

number, since y r w / p is a velocity (called the friction velocity). The solution of equations (4-75) and (4-76), 
with boundary conditions at x 2 = 0 and 2(x 2 ) c , gives 
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^ + ,v + 

u. ~ u. 

1 1 


+ / +\ 
X k ’ ’ '2 ' c 


and 


+ + 
a = a 


+ ,+ / +\ 

x k , t , (x 2 ) 


2 'c 


However, since mean values are steady and fully developed, u*u* (see eq, (4-14)), and u* = U + (see 
eq. (4-18)), are not functions of x|, Xg, or t + . Thus, 


and 


U 1 U 2 = U 1 U 2 


’ ( X 2)c 



(4-78) 


(4-79) 


Next let (x^) c — > «>. Then for the region where x 2 << (x 2 ) c , t * ie parameter (x 2 ) c drops out of 
the determining boundary conditions. That occurs because the only boundary condition determining the 
turbulence at x 2 << (x 2 ) c is the one at the near wall, the boundary condition at x 2 = 2(x 2 ) c — > «> 
having no effect. So equations (4-78) and (4-79) become, for (x 2 ) c —> <» and x 2 « ( x 2 ) c > 


and 


+ + + + , +. 
U 1 U 2 = U 1 U 2 ( X 2^ 


(4-80) 



(4-81) 


For the present fully developed unidirectional flow, the eddy viscosity from equation (4-28) is 


e 


u,u 


1 u 2 


dU 1 /dx 2 


or 


So, from equations (4-78) and (4-79), 



+ + 
U 1 U 2 


dl^/dx* 


S S + , +, 

— — — X 2 ’ ( X 2 )c > 

V l / 1 J 


(4-82) 


(4-83) 


(4-84) 
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or, from equations (4-80) and (4-81), 



(4-85) 


for x+ « (x+) c and (x+) c co. 

The fact that equations (4-80), (4-81), and (4-85) hold for x 2 « (x 2 ) c and (x 2 ) c — ♦ «> is often 
called the law of the wall. Although that is usually considered to be an empirical law, we have obtained 
it here by using the Navier-Stokes equations in a thought experiment. Note that if one considers the 

region so close to the wall that u*u* — > 0, the so-called laminar sublayer, equations (4-71) and (4-77) 
give (for x 2 « (x^)J 


or 




(4-86) 


(4-87) 


4.3. 2. 7 The velocity-defect law.— Consider now the central region of the channel, sometimes 
referred to as the region away from the wall. As in the last section the flow is fully developed and 
statistically steady. We write equations (4-75) and (4-76) as 


and 


X _ ^ d(ujX) 

dt* dx? 

k 


d<r i ^ u i 

X ( X 2 ) c KK 



dx*dx* dx*dx? 

it i k 


where 


(4-88) 


(4-89) 



(4-90) 


The quantities u| and cr + are defined in equation (4-77). Proceeding as for the law of the wall in the 
last section, we note that equations (4-88) and (4-89) constitute a set of four equations in the four 
unknowns Uj , u 2 , U 3 , and a + . Here the walls where the boundary conditions are set are at x 2 = 0, 2 
(e.g., fit = 0 at x 2 = 0, 2 ). As in the last section we do a thought experiment, starting with given 
initial conditions. We consider a cross section where the mean flow is fully developed and continue the 
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solution until the flow is statistically steady. Then equations (4-88) and (4-89), or their solution, show 
that 


and 


„ + „ + 
u. = u. 

1 1 


*;> (x 2 ) c 


* + ~ + 
a — a 




(4-91) 


(4-92) 


Note that whereas in the last section (x 2 ) c was included in the functional relations because of its use for 
the location of the boundary condition at x 2 = 0, 2(x 2 ) c , here it is included because of its appearance in 
the differential equation (4-88); (x 2 ) c is not needed for the location of the boundary condition at x 2 — 2. 

* * ^* nce * mean values are steady and fully developed, u 2 = U 1 (see eq. (4-18)) is not a function of 

x l> x 3 ? or t , as are the unaveraged quantities in equations (4-91) and (4-92). Thus 




(4-93) 


Now we confine our attention to the region away from the wall and let (x 2 ) — > «>, so that the vis- 
cous term, including (x 2 ) c , drops out of equation (4-88). Batchelor explains this loss of the influence of 
viscosity in the context of Fourier analysis (ref. 12). He notes that the region of wavenumber space 
which is affected by the action of viscous forces moves out from the origin toward a wavenumber of infin- 
ity as the Reynolds number increases. In the limit of infinite Reynolds number (infinite (x 2 ) c in our 
case) the sink of energy is displaced to infinity, and the influence of viscous forces is negligible for wave- 
numbers (reciprocal eddy sizes) of finite magnitude. (See chapter V for a discussion of Fourier analysis.) 
Thus, for (x 2 ) c — > co, equation (4-93) can be written, for the region away from the wall, as 


w - «)c = - K), 




(4-94) 


where (U^) c is the value of U| at the channel midpoint. Equation (4-94) is written with LTj — (U)() c 
as the dependent variable to ensure that curves for different Reynolds numbers will collapse at the chan- 
nel midpoint. According to equation (4-94) they will then also collapse at points near the midpoint. The 
quantity U| — (U^) c is called the dimensionless velocity defect, and equation (4-94), which applies in 
the central region of the channel when viscous forces are negligible there (when (x 2 ) c — > co), is known as 
the velocity-defect law. 1 As in the case of the law of the wall, equation (4-94) is usually considered to be 
an empirical law, but it is obtained here by using the Navier-Stokes equations in a thought experiment. 


1 Note that letting (x 2 ) c — > co and neglecting the viscous terms in equations (4-88) destroys the 
accuracy of the solution near the wall, since viscous effects are not negligible there. But as already 
shown, letting (x 2 ) c — > co preserves the accuracy in the region away from the wall, the region we are 
interested in here, since viscous effects are negligible away from the wall. 
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4. 3. 2*8 The logarithmic law . — Next we want to determine the form of the mean velocity profile 
required in a possible overlap region, where equation (4-81) (the law of the wall), and equation (4-94) 
(the velocity defect law), both apply. From equations (4-71), (4-83), and (4-84), which apply over a 
whole channel cross section in the fully developed region, we get 


+ 






dx 2 


(4-95) 


It is known that, except in the immediate vicinity of the wall, where x 2 is less than say 70, the 
molecular shear stress is much smaller than the turbulent. We will show later in this section that onr 
results are consistent with that statement. Thus, neglecting the molecular shear stress in equation (4-95) 
we have 


+ 






(4-96) 


In order that equation (4-81) (the law of the wall) will hold, we confine the range of equa- 
tion (4-96) to x 2 /(x 2 ) c << 1, but at the same time keep x 2 /(x 2 ) c large enough that the turbulent 
stress is large compared with the molecular stress. 2 For that range we can use equation (4-85) in equa- 
tion (4-96). The latter then becomes 


1 


-W 


dx„ 


(4-97) 


Note that equation (4-97), when integrated from x 2 = 0 (at the near wall) to x 2 , is in the form of equa- 
tion (4-81), and is thus an expression of the law of the wall. 


To ensure that equation (4-97) will also satisfy equation (4-94), the velocity-defect law, we first 
write the former in terms of x * = x 2 /(x 2 ) c = x 2 /(x 2 ) c , (x 2 ) c5 and the dimensionless velocity defect 
(Uj) c — Uj\ That gives 



e 

i/ L 


/ +\ * 
( X 2 )c X 2 


.«) c - u; 

(x 2 ) c dx* 



(4-98) 


2 The range of x 2 for which that is possible increases as (u 2 ) c 
obtaining the velocity-defect law and the law of the wall. 


oo; the latter was assumed for 
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where the subscript c refers to values at the channel center. If equation (4-98) is to agree with equa- 
tion (4-94), then (x 2 ) c must cancel out of equation (4-98). For that to happen, sjv must be propor- 
tional to (x 2 ) c , so that 


e , 

-(xj 


e 

v L 


( X 2 )c x 2 


= K ( x 2 )c x 2 = Kx * ’ 


(4-99) 


where K is a constant (called the Karman constant). Then equation (4-98) becomes 



(4-100) 


• Sfc £ 

which, when integrated from x 2 = 1 (at the channel center) to x 2 , is in agreement with equation (4-94). 
Thus, only when sjv is given by equation (4-99) does equation (4-97) satisfy both equations (4-81) 
and (4-94), as it must in an overlap region. 3 Substituting equation (4-99) into equation (4-97) gives 


or, on integration, 




> 


U. + = 2llnx* + C . 
1 K 2 


(4-101) 


(4-102) 


Equation (4-102) is the well known logarithmic velocity distribution which applies to the portion of the 
law-of-the-wall region where the turbulent shear stress is much greater than the molecular. It can be 
obtained in a number of ways (e.g., ref. 15), but is obtained here from equations (4-81) and (4-94) by 
assuming an overlap region. That was first done by Millikan, but in a somewhat different way (ref. 16). 
Besides obtaining the logarithmic law, we have shown that any expression for e/v which satisfies both 
the law of the wall and the velocity defect law in some overlap region must, according to our analysis, 
reduce to Kx 2 in that region (eq. (4-99)). 

We can now see whether the assumption that the turbulent shear stress is much greater than the 
molecular shear stress except very close to the wall, where x 2 is less than about 70, is in agreement with 
our results. Substituting the expression for sjv in equation (4-99) into equation (4-95) gives 


3 

Actually, there are expressions for sjv which satisfy both equations (4-81) and (4-94) but which 
have an appearance which differs from that of equation (4-99). However, they reduce to equation (4-99). 
For example, Prandtl’s mixing-length expression (ref. 6) and von Karman ’s similarity expression (ref. 15) 
are both of the form ( sjv ) = (e/i/)(U|, x 2 ) and satisfy equations (4-81) and (4-94). But (e/v),= 

(*/*') (^l> x 2) = (*/*') [Uj (x+), x 2 ] = (e/i/)(x 2 ). This agrees with equation (4-85) and leads to equation 
(4-99). Equation (4-99) is simpler than the other expressions for e/v and, since the other expressions 
reduce to it, may be more fundamental. 
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(4-103) 
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+ 



= (1 + Kx 2 + ) 



so that the ratio of turbulent to molecular shear stress is Kx^/1. Since K « 0.4 the turbulent shear 
stress is more than an order of magnitude greater than the molecular stress for x 2 > 70. 


4.3.2.9 An expression for e/v for fully developed flow which applies at all distances from a 
wall . — We now give a general expression for the eddy diffusivity which satisfies conditions given in 
sections 4.3.2.4 and 4.3.2.6 to 4.3.2.8, and which applies across the whole channel in the fully developed 
region. The expression is 


where 


e r + 

. = Lx, 


( x 2 *ic 


5/4 


(4-104) 


L = 


+ 13/5 13/6 

allj for aU 1 < 2/5 

, 13/6 

2/5 for aU* >2/5 


(4-105) 


and a = 3x10 4 . 

Using equation (4-104) to close equation (4-95), a fully developed form of the equation for the mean flow, 
we get 


+ 
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5/4 

1 + Lx 2 + 

1 - 2 



(*2 + )c 



du; 


dx„ 


(4-106) 


where L is given by equation (4-105). 

Consider now whether the conditions obtained in sections 4.3.2.4 and 4.3.2.6 to 4.3.2.8 are satisfied 
by equations (4-104) to (4-106). First note that those equations satisfy equations (4-79) and (4-84) for 
0 < x+ < oo. Then for *2 « K) c they agree with the law-of-the-wall equations (4-81) and (4-85). As 
x 2 0? equation (4-86) or (4-87) is satisfied. That is, as x 2 — ► 0, Uj becomes equal to x 2 . An addi- 
tional condition satisfied as x 2 — > 0 is that e/v — > x 2 * 8 /^, which is equivalent to equation (4-51). This 
follows from evaluating equations (4-104) and (4-105) for small values of x 2 and using equation (4-87). 
Then by using equations (4-83) and (4-86), we get equation (4-51), which in turn satisfies equations (4-46) 
to (4-50). These last equations are consequences of continuity and the nonslip boundary condition at the 
wall, so the expression for e/v given by equations (4-104) and (4-105) is also in agreement with those 
consequences. 
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Next consider whether equation (4-106) satisfies the velocity defect law (eq. (4-94)) for the region 
away from the wall, where the turbulent shear stress is much greater than the molecular (where x 2 is 
greater than say 70). For that region equation (4-106) can be written as 


2 * 
X 2 


X„ (1 


- O' 1 * 


(«D. - u i* 


dx„ 


(4-107) 


where x 2 = X 2 /(x 2 ) c . But equation (4-107) when integrated from the channel center to x 2 , is in the 
form of equation (4-94). Thus equation (4-106), from which equation (4-107) is obtained for the region 
away from the wall, satisfies the velocity-defect law in that region. Integration of equation (4-107) gives 


K) c 
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1 + (1 — x 2 ’)‘/ 4 
i - (i - x ;) 1 / 4 
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(4-108) 


for our velocity-defect law. 

Finally, for the law-of-the-wall region, where X 2 « (x^) c , equation (4-106) becomes 


1 = 


i + * x ; 

5 2 


du; 


(4-109) 


dx„ 


Then for the part of that region where the turbulent shear stress is much greater than the molecular (the 
region away from the wall), we have 


which, when integrated, becomes 



IL + = 1 In x„ + + C , 
1 2 1 


(4-110) 


(4-111) 


which agrees with the logarithmic distribution in equation (4-102). 

Equation (4-106) is plotted semilogarithmically for several values of (x 2 ) c in figure 4-2. Also 
included in the plot are some experimental data from references 11 and 15, and a numerical solution 
(direct numerical simulation) of the unaveraged Navier-Stokes equations from reference 17. Although the 
experimental data are for a pipe, and equation (4-106) was obtained for a flat channel, a similar equation 
is obtained for a pipe by letting (x 2 ) c be the radius of the pipe. For the region close to a wall, the flow 
in a pipe should be similar to that in a channel, since, for that region the pipe wall could be considered 
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flat. Also, for a pipe, as well as for a channel, the total shear stress varies linearly with wall distance, as 
can be shown by writing a force balance on an element of fluid for both configurations. It is conceivable 
that the eddy diffusivity for the central region might be different for the two configurations, but such 
differences are not apparent in figure 4-2. Equation (4-106) appears to represent both the data and the 
numerical solution quite well. 

We can calculate friction factor as a function of Reynolds number by using the velocity distribu- 
tion from equation (4-106). Since we will be comparing the calculated friction factors with those obtained 
experimentally for a pipe, the quantity (x 2 ) c will be the pipe radius (see the last paragraph). The fric- 
tion factor f is given by 

pu> u ; 2 

and the Reynolds number Re is obtained from 

Re = 2 (x 2 ) c ^ = 2 (x 2 *) c U;, (4-113) 


where U a is a dimensionless bulk or mixed-mean velocity and is given by 
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(4-114) 


A plot of friction factor versus Reynolds number, as obtained from equations (4-106), (4-112), 
and (4-13), is compared with experimental data for turbulent flow in a pipe in figure 4-3.' As in the case 
of the velocity distribution in figure 4-2 the agreement is quite good. 


4.3.2.10 Thermal law of the wall, temperature-defect law, and temperature logarithmic law. — Laws 
analogous to those obtained for the mean velocity distribution in sections 4.3. 2.6 to 4.3. 2. 8 can be 
obtained for the mean temperature distribution when fully developed heat transfer, in addition to fully 
developed flow, occurs in a passage. The procedures are similar, and so this section will not be detailed. 
The main difference here is that the energy equation (4-12) (in appropriate dimensionless forms) is added 
to the lists of unaveraged equations used in the thought experiments. Thus, in addition to equa- 
tions (4-75) and (4-76) we use 


dT + a (T + u k ) i 
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and in addition to equations (4-88) and (4-89) we have 
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~ + _ (T w - T)cr w v 

where T = , Pr = _ is the Prandtl number, and (Q 2 ) w is the time-averaged normal 

(Q2)w \j T v// P 

heat transfer at the wall. The rest of the dimensionless quantities have already been defined in equa- 
tions (4-77) and (4-90). Then one obtains the thermal law of the wall, the temperature-defect law, and 
the temperature logarithmic law in place of equations (4-81), (4-94), and (4-102) respectively. These are 

T + = T + (x 2 + , Pr), 


(T + - T + ) = (T + - T + )(x!), 


and for an overlap region, where both the thermal law of the wall and the temperature defect law apply, 

T + = J_ In y + + C(Pr) , 

K l 


where T + = T + is the time-averaged dimensionless temperature difference, and Kj = K (see 

eq. (4-102)) if = e as in equation (4-38). Note that the Prandtl number Pr does not appear in the 

temperature-defect law because the term 1 / [(x^ )Pr]| d 2 T + / 3x*dx* drops out of the dimensionless 

energy equation used for the temperature-defect law when (xj) — > oo. (This loss of the influence of ther- 
mal conduction (or of thermal smearing) can be explained by an argument similar to that for the influ- 
ence of viscosity, as given after equation (4-93).) On the other hand, Pr appears in the thermal law of 

the wall because the term (l/Pr)d 2 T + /dx. + dx. + does not drop out of the dimensionless energy equation 

used for the law of the wall when (xj),. — > °°. Finally, note that in getting the logarithmic law for the 
temperature distribution we use, as obtained from equations (4-29) and (4-74) 

q 2 


(Q 2 )v 


— + — 
Pr v 


dT 


(4-115) 


dx„ 


in place of equation (4-95), where is the eddy diffusivity for heat transfer. As a result of this change 
s is replaced by e h and/or Uj by T + in equations (4-97) to (4-102). Note that, as for equation (4-95), 
and except for very small Prandtl numbers, if one lets (x^),, — > 00 there will be an extensive region where 
the turbulent term is much greater than the molecular term, and where the term on the left side can be 
replaced by 1. The results in this section, like those in sections 4. 3. 2.6 to 4. 3. 2. 8 follow from a thought 
experiment with the unaveraged continuum equations. 

4.2.3.11 A calculation of fully developed convective heat transfer . — Consider next a treatment of 
the fully developed heat transfer which can occur when temperature gradients exist in a flowing fluid. If, 
in accordance with equation (4-38) we set = e, equation (4-115) becomes 
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where, again, 


and 


Q2 

(Q 2 )w 
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(Q2)w ft . 



(4-116) 


is the Prandti number of the fluid. As for the velocity distribution (eq. (4-106)), e/v in equation (4-116) 
is obtained from equation (4-104). Equation (4-116) applies to either a channel or a pipe, as does equa- 
tion (4-106). However, the quantity Q 2 /(Q 2 ) W in equation (4-116) is different for the two configura- 
tions, although the effect of the difference on the heat-transfer coefficients is small except for liquid 
metals (which are not considered here). Since we will be comparing our results mainly with experimental 
heat-transfer (and mass-transfer) data for a pipe, we use Q 2 /(Q 2 ) W for a pipe (ref. 7 or 9): 


q 2 

(Q 2 )w 










(4-117) 


The quantity (x 2 ) c is defined as (x 2 ) c y r u /p /V, where (x 2 ) c is now the pipe radius, U| is obtained 

from equation (4-106) and U+ is given by equation (4-114). (Compare eq. (4-117) with that obtained 
for a channel by using the first and last members of eq. (4-74)). 


We want to calculate the Stanton number as a function of Reynolds and Prandti number. The 
Stanton number is given by 

h 1 
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where the heat-transfer coefficient h is defined as 
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T - T 
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and the dimensionless mixed-mean temperature difference T* is obtained from 
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where U ^ is calculated from equation (4-106), and T + is obtained from equations (4-116) and (4-104). 
Finally the Reynolds number Re is given by equation (4-113). By varying the parameter (x 2 ) c , Stanton 
number can be calculated as a function of Reynolds number and Prandtl number from these equations. 

Figure 4-4 shows plots of calculated Stanton number versus Prandtl number for three Reynolds 
numbers. Included for comparison are experimental data for heat and mass transfer. The mass-transfer 
data are included by replacing temperatures by concentrations and molecular thermal diffusivities by 
molecular diffusivities for mass transfer. The agreement of the experimental data with the solid curves 
calculated by using velocity and temperature distributions from equations (4-106), (4-116), and (4-104) is 
good over the entire range of variables shown. 

Of some interest is the fact that the good agreement with experiment was obtained by setting the 
eddy diffusivity for heat transfer (or for mass transfer) equal to that for momentum transfer e, as 
obtained from the mixing-length theory (see eq. (4-38)). By contrast, the agreement obtained by setting 
e / e h ec l ua l to the molecular Prandtl number (about 0.7 for air), as is often done, is much poorer. This is 
indicated by the dashed curve in figure 4-4c. These results do not, of course, prove that e/e-^— 1 at all 
points in the flow, since e /e^ may vary with position. They only show that the mean effective value of 
e/ejj is very close to 1. 

4.3.2.12 Some other closure assumptions for fully developed or nearly fully developed (equilibrium) 
flows . — The results in the last section utilized an approximate expression for the eddy viscosity which 
was, however, based partially on exact information. Here we review briefly some earlier approaches. 

A reasonable expression for the region away from walls is the von Karman similarity expression 
(ref. 33). That expression is most easily obtained by assuming that, away from boundaries, the turbu- 
lence at a point is a function only of conditions in the vicinity of the point, in particular, of the first and 
second derivatives at the point. Then by dimensional analysis we obtain for the eddy viscosity 
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(4-121) 


where K is the Karman constant. 

Von Karman’s hypothesis (eq. (4-121)) has been proposed for the region away from a wall. Close 
to a wall we assume that the eddy viscosity is a function only of quantities measured relative to a wall, 

Uj and x 2 , and of v. The simplest assumption consistent with dimensional analysis and the require- 
ment that the effect of v should become small for large x 2 is then (ref. 33) 

(4-122) 
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where n is an experimental constant (n = 0.124). Equations (4-121) and (4-122) give results for flow 
and heat (or mass) transfer in tubes which are almost as good as those in figure 4-4 (see ref. 33). The use 
of those equations also gives good results for the boundary layer on a flat plate (ref. 34). 


By making a small modification in the derivation of equation (4-121) we can obtain an expression 
for e which should be applicable to a vortex. We assume that the turbulence at a point is dependent 
only on the shearing deformation rate at the point and in the vicinity of the point-, that is, it is a function 
of the deformation rate and its derivatives. If we exclude derivatives of the deformation rate higher than 
the first we get, for a parallel flow (eq. (4-121)). However, for a circular vortex, we get 
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(4-123) 


where v is the tangential velocity and r is the radius. For a vortex flow, v = v Q r 0 /r, where the sub- 
scripts 0 refer to values at some arbitrary radius. Equation (4-123) then becomes 


K 2 v 0 r 0 

s = 

2 


(4-124) 


Equation (4-124) has been used profitably for turbulent vortex flows 
physical clouds (ref. 36), and in the atmosphere (ref. 37). 

Finally we mention closures based on Prandtl’s mixing-length 

" 1*2 = -«2 



in vortex tubes (ref. 35), in astro- 
? 2 (ref. 6), where { 2 is defined by 

(4-125) 


Prandtl assumed that ? 2 = Kx 2 , where x 2 is the distance from a wall. Van Driest (ref. 38) has modi- 
fied that assumption by introducing a damping factor which reduces J 2 i n the region close to a wall: 


< 2 = Kx 2 

where A is an additional experimental constant. Equation (4-126) appears to be reasonably applicable 
to the regions both close to and away from a wall, as is our equation (4-104). Equation (4-104) gives, in 
addition, a velocity profile which is accurate in the so-called wake region, where x 2 ~ (x 2 ) c (see fig. 4-2). 
Although the use of an equation which is applicable to two or more regions in a flow may be somewhat 
more convenient than the use of a separate equation for each region, there is basically no reason why one 
equation should apply to more than one region; the turbulence mechanism is likely different in the dif- 
ferent regions. 


1 — exp 


(4-126) 
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A word should perhaps be said about the value of the Karman constant K. If ut = uJK, 

(x 2 )], as in figure 4-2, then a value of K = 0.4, or 2/5, (see eq. (4-105)) gives good results. On the 
other hand, if it is assumed that U| = U{(x 2 ), so that a single curve is obtained instead of the multiple 
curves in figure 4-2, then K = 0.36 gives better agreement with the data (ref. 33). Finally, for a vortex 
(eq. (4-124)), K seems to be close to 0.3 (ref. 35). It is not surprising that the value of K appears to be 
slightly different for a vortex than for a parallel flow, since, as for all closures, the ones considered here 
are approximate and correspond only partially to reality. Thus the constants in a closure scheme must 
often be “fine-tuned” when the flow configuration is changed. 

4.3.2.13 A treatment of moderately short highly accelerated turbulent boundary layers. — The last 
few sections considered closure schemes for fully developed or nearly fully developed turbulent flows. We 
turn now to a problem which is in some ways the opposite of those just considered: Given a fully devel- 
oped or nearly fully developed turbulent flow and transverse heat transfer to which the boundary-layer 
assumptions are applicable, to determine the streamwise evolution when the flow is subjected to a severe 
streamwise pressure gradient. A boundary layer subjected to a severe pressure gradient can be considered 
to be the same thing as a highly accelerated boundary layer. According to Bernoulli’s equation, which 
applies outside the boundary layer, where viscous and turbulence stresses are negligible, 
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where is the velocity just outside the boundary layer. For a thin steady-state two-dimensional 
boundary layer with constant properties and without buoyancy, equations (4-25), (4-26), and (4-21) 
become respectively 
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(4-129) 


(4-130) 


where Xj and x 2 are respectively in the direction along and normal to the wall. The pressure gradient 
is written as a total derivative dP/dx 1 because P is not a function of x 2 for a thin boundary layer. 
Equations (4-128) to (4-130) apply to a boundary layer even when the flow is along a curved wall 
(ref. 39). 

In order to solve equations (4-128) to (4-130) to obtain the evolution of and T, u^u 2 and 
ru 2 must be known at each point in the flow. The experiments of Blackwelder and Kovasznay (ref. 40) 
suggest that although severe pressure gradients caused the mean flow in those experiments to change con- 
siderably along streamlines, the Reynolds stresses, at least in the important intermediate region of wall 
distances, were relatively unaffected. This leads us to the hypothesis of frozen Reynolds stresses and 
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turbulent heat transfer along streamlines in a moderately short highly accelerated turbulent boundary 
layer. We gave theoretical arguments for that hypothesis in references 41 and 42, where we showed that 
it gives good agreement with experiment when applied to equations (4-128) to (4-130). Launder (ref. 43) 
had earlier used the concept of a frozen Reynolds stress, but in an approximate integral approach to the 
boundary layer rather than in the solution of the partial differential equations (4-128) and (4-130). 

Changes along streamlines, if not zero, should be smaller than say changes at constant x 2 , since flow 
along streamlines tends to involve the same fluid or eddies at various streamwise points. In applying the 
hypothesis of frozen UjU 2 and ru 2 along streamlines it is convenient to transform equations (4-128) 
to (4-130) from (x 1? x 2 ) to (x^ ip) coordinates (von Mises coordinates), where the stream function ip is 
given by 
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where the subscripts are quantities held constant. Equations (4-128) to )4-130) then become 
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where the subscripts ip indicate changes along a streamline (at constant ip)- By virtue of equa- 
tion (4-131), equations (4-133) and (4-134) satisfy continuity automatically. 


The equations for the evolution of UjUj and of ru 2 will be obtained later in this chapter (see 
eqs. (4-140) and (4-141)). For our purpose the important point is that those equations do not contain 
dP/dxj, in contrast to equation (4-133) for the evolution of Uj or to equation (4-127) for the evolution 
of U^, the latter of which is a special case of equation (4-133). Taking that fact into account, as well as 
the discussion in the paragraph preceding equation (4-131), one would expect, for large dP/dx 1? that there 
would be large gradients of U w or of Uj along streamlines without correspondingly large streamwise 
gradients of UjU 2 and ru 2 . Then if the boundary layer is moderately short, the hypothesis that the 
turbulent shear stress and transverse heat transfer are frozen at their initial values as one proceeds along 
a streamline in a flow with large streamwise gradients of or P should be a good approximation. 

Rescaling the variables in equations (4-133) and (4-134) so as to convert them to dimensionless 
form, and introducing the frozen turbulent stress/heat-transfer hypothesis give 
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and the arrow — > means “has been replaced by.” The quantity x 0 is the value of x^ at the initial 
station, U Q is the velocity outside the boundary layer at the initial station, (T w ) 0 is the wall tempera- 
ture at the initial station, is the constant temperature outside the thermal boundary layer, and Pr 
is the Prandtl number. The quantities (u 1 u 2 ) 0 and (ni 2 ) 0 are, respectively, the values of u^u^ and 
ru 2 on the same streamline as u 1 u 2 and ru 2 but at the initial station. 

The frozen turbulent stress/heat-transfer hypothesis can be tested by numerically integrating equa- 
tions (4-135) and (4-136) along streamlines (at constant if)) and comparing the results with experiments 
for the same conditions. Upstream of the regions of severe streamwise pressure gradients the pressure 
gradients were essentially zero. 

Velocity profiles (Uj/U^ against x/j/u) are plotted and compared with experiment in figure 4-5. 
(Note the shifted vertical scales.) In all cases the effect of the pressure gradient and the total normal 
strain parameter U^/Uq is to flatten the profiles. The agreement between theory and experiment is 
considered good. 

Semi-logarithmic plots of Uj/(r w /p) 1 / 2 against (r w /p) 1 / 2 x 2 / v (law-of-the-wall plots) are given in 
figure 4-6. These profiles show the inner region of the boundary layer much better than does figure 4-5. 
The shear stress r w at the wall for the theoretical curves was obtained from the slope of the velocity 
profile at the wall by using points very close to the wall ((~ w //>) 1//2 x 2 /v << 1). Points very close to the 
wall were necessary because of the nonlinearity of the profile close to the wall in the presence of a severe 
pressure gradient. It might be pointed out that this nonlinearity makes the experimental determination 
of the shear stress at the wall extremely difficult. Both theory and experiment indicate that the original 
logarithmic and wake regions are destroyed by the pressure-gradient and normal-strain effects, although a 
new logarithmic layer of smaller slope seems to form eventually. Also, the thickness of the sublayer 
approximately doubles, indicating an apparent “relaminarization,” as observed experimentally by many 
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investigators. However, it is not a true relaminarization since, at least in the theory, u,u 2 is constant 
along streamlines. The agreement between theory and experiment appears to be quite good. For values 
°f UJU 0 larger than those shown, the approximation of a frozen Reynolds stress apparently begins to 
break down. 


In order to see how sensitive the development of the mean profile is to the Reynolds shear stress, 
results were calculated for UjU 2 = 0 and are shown dashed in figures 4-5(b) and 4-7. The effect of u^u 2 
on the profiles in figure 4-5(b) is slight. The law-of-the-wall plots in figure 4-7, on the other hand, show 
a significant quantitative effect of UjU 2 on the profiles, but qualitatively the curves for UjU 2 == 0 and 

u l u 2 * 0 are muc h the same. In both cases the original logarithmic and wake regions are destroyed and 
the sublayer is thickened. The difference between the indicated quantitative effects of ^u 2 on the pro- 
files in figures 4-5 (b) and 4-7 is evidently due to the difference in scales and in scaling parameters in the 
two figures. 

Figure 4-8 shows, for a large value of the pressure-gradient parameter, the contributions of 
various terms in equation (4-128) to the rate of change of the nondimensional mean kinetic energy 
(l/2)d(Uj/U 0 ) / ^(x^Uq/ i/). (The energy in the transverse velocity component is negligible for a bound- 
ary layer.) The contribution of the Reynolds stress term is very large in a narrow region near the wall. 

However, that tends to be offset by the viscous contribution. Comparison of the curves for UjU 2 — 0 

and u 1 u 2 * 0 shows that the viscous contribution adjusts its value so as to offset the effect of the 
Reynolds stress. The viscous term is not zero at the wall but balances the pressure-gradient term so that 
Uj can remain zero at the wall. Thus, the present velocity profile, in contrast to the case of zero pres- 
sure gradient, is nonlinear at the wall. (If it were linear the viscous term in equation (4-128) would be 
zero at the wall.) 

The pressure-gradient term is independent of wall distance and, for the case shown in figure 4-8, 
becomes dominant for (r w /p)^'^x 2 / v > 40. Thus, the destruction of the logarithmic and wake regions is 
due mainly to the pressure-gradient term, rather than to a change in the structure of the turbulence 
(although some change in structure may occur (ref. 45)). Also, the thickening of the sublayer is mostly, 
although not entirely, due to the pressure-gradient contribution, since as mentioned, viscous effects tend 
to offset the Reynolds stress contribution. 


The results in figure 4-8 are, of course, for a large pressure-gradient parameter. For regions of 
lower pressure gradient, the Reynolds stress will have a greater effect, as shown in figure 4-7. Also, the 
velocity profile at any position depends on the whole distribution of pressure gradients up to that posi- 
tion; that is, Uj is a functional of dP/dxj, or 


U, = U, 


dP 

dx. 


(€) 



(4-138) 


where 0 < £ < Thus, there is a quantitative (but not a qualitative) effect of UjU 2 on the velocity 
profile, even at those positions where the pressure-gradient parameter is large. 

The analysis can be easily extended to include mass injection at the wall by transforming equa- 
tion (4-135) from (xj, VO to (xj, i/t ') coordinates, where ip 1 = if> — V> w , and V w * s the stream function 
at the wall. The latter will vary with Xj in accordance with equation (4-131). Equation (4-135), when 
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written in (xj, ip ) coordinates, has the additional term — (l/2)(U 2 ) w 6Uj/ dip' on the right side, where 
( U 2) w is the dimensionless normal velocity at the wall, and where, as in equation (4-135), all quantities 
^ ave k een nondimensionalized by Uq and v (see eq. (4-137)). As before, we use the simplification that 

u l u 2 remains frozen as it is convected along streamlines (along lines of constant ip, not constant ip'). 
The injected fluid is assumed to be turbulence-free. 

To show the effect of mass injection on a boundary layer with severe pressure gradients, mass 
injection was added in the theoretical calculations in figure 4-6. For positive mass injection figure 4-9 
shows that the normal flow quickly raises the U J{tJp) 1 I 2 curve, particularly in the wake region, after 
which the favorable pressure gradient lowers and flattens the curve. The resulting curve still lies above 
the initial profile. For negative injection, the normal flow and pressure gradient lower and flatten the 
initial profile. These trends are similar to those observed in the experiments of Julien, Kays, and Moffat 
(ref. 46). 


All of the results so far were for favorable pressure gradients, but the analysis should apply as well 
to severe unfavorable gradients. Figure 4-10 shows a comparison between theory and experiment for the 
results of Kline et al. (ref. 45) for a severe unfavorable pressure gradient (their fig. 9(a)). The results 
indicate that the adverse pressure gradients produce an exaggerated wake region, but that the logarithmic 
and sublayer regions are relatively unaffected. The agreement between theory and experiment is good. It 
might be mentioned that the results for adverse pressure gradients were more sensitive to the UjU 2 dis- 
tribution than were those for favorable gradients. In particular, when UjU 2 was taken as zero, sepa- 
ration occurred upstream for the run shown in figure 4-10 for U^/Uq = 0-92. Thus the presence of 
turbulence appears to delay separation. This is evidently because the Reynolds stress term in equa- 
tion (4-128) is positive close to the wall and thus tends to increase U x in that region. 

Figure 4-11 shows theoretical Stanton-number and skin-friction-coefficient variations with dimen- 
sionless longitudinal distance for run 12 from reference 47. The shear stress and heat transfer at the wall 
for the theoretical curves were obtained from the slopes of the velocity and temperature profiles at the 

wall by using points very close to the wall y (r w / p )x 2 / ^■Cl) . Also included in the plot are experimental 
values for U^/Uq- Initial conditions were taken at x 0 = 4'32 ft. The difference between the Stanton 
number and skin-friction-coefficient variations is rather striking and indicates that Reynolds analogy (see, 
e.g., ref. 33) does not apply in regions of severe pressure gradients. This difference is also indicated in the 
experimental results of references 40 and 47. 

The effect of favorable pressure gradients on velocity and temperature distributions is illustrated in 
figure 4-12, where theoretical values of T + and are plotted against x 2 for a low and a high value 
of pressure-gradient parameter K = {yj pU^ ) )dP/dx 1 . The results are again for run 12 from reference 47. 
The effect of the pressure gradient on the T + profile tends to be opposite to that on the U 2 profile. 
Whereas the pressure gradient flattens the profile in the outer region of the boundary layer, it 
steepens the T + profile in that region. The same trends have been observed experimentally in refer- 
ences 44 and 48. The difference between the velocity and temperature results (or the skin friction and 
heat transfer results) is evidently due to the fact that the equation for the evolution of the mean velocity 
contains a pressure-gradient term (eq. (4-128)) whereas that for the mean temperature does not (eq. (129)). 
Although the temperature equation does not contain a pressure gradient term, the pressure gradient can 
still affect the temperature through the mean velocity. 

A comparison between theory and experiment for the evolution of Stanton number in severe 
favorable pressure gradients is presented in figure 4-13 for three values of maximum pressure-gradient 
parameter K m . The value of x Q (the initial station for each run was taken at the point where the local 


37 



pressure-gradient parameter K L starts to increase rapidly). For the smaller values of K m good agree- 
ment between theory and experiment is indicated for values of U^/Uq which are not too large. It 
appears that the range of values of U^Uq for which the theory applies increases as the pressure- 
gradient or acceleration parameter increases. For values of U Q (xj - x Q )/i/ (or of U ^/Ug) greater than 
those shown in figure 4-13, the theory appears to break down because the total streamwise strain becomes 
too great and/or the local pressure-gradient parameter K L becomes too small. 

Figure 4-13 also shows Stanton numbers calculated for turbulent initial velocity and temperature 
profiles, but for ru 2 = u x u 2 = 0 . It is seen that the turbulent stresses and fluxes have a very large 
effect on the evolution of the Stanton number. The large effect of turbulence on Stanton number perhaps 
raises questions about calling regions such as those shown in figure 4-13 “relaminarization regions,” 
although at the large values of K m and U Q (x 1 — x^)/v there is some tendency for the zero-turbulence 
curves to approach those for turbulence. The effects of turbulence on the velocity and temperature dis- 
tributions are also considerable, although the general trends without turbulence are similar to those in 
figure 4-12. 

It should, of course, be remembered that, as is the case for Uj (eq. (4-138)), Stanton number, as 
well as C f , and T at any longitudinal position depend on the whole distribution of pressure gradients up 
to that position. That is, Stanton number, for instance, is a functional of dP/dx 1? or 

(4-139) 

where 0 < £ x r 


St = St 


dP 

dxi 


(0 


4.3.3 One-Point Correlation Equations 

Sections 4.3.1 and 4.3.2 considered the equation for the evolution of the mean flow (eq. (4-25)) and 
that for the mean temperature (eq. (4-26)). Those were obtained by averaging the equations for the 
instantaneous velocity and temperature (eqs. (4-11) and (4-12)) after applying Reynolds decomposition, 
and contained the important but undetermined turbulent stress tensor — pUjUj and turbulent heat trans- 
fer vector pcujf. This was our first encounter with the closure problem, and sections 4.3.2 to 4.3.2.13 
considered some simple closure schemes. 

One can construct equations for the evolution of the undetermined quantities UjUj in equations (4-25) 
and (4-27) from equation (4-22) for gj = 0 and the following similar equation for the component Uj: 


dUj 

dt 


dx 


d , s 1 do ^ ^ u j 
(UjU k ) + v. 


- Ut, 
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u k— + 
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UjU k 


Note once again that <r g drops out the equations of motion for gj = 0 by virtue of the equation follow- 
ing equation (3-22). Next multiply equation (4-22) by Uj and the above equation for Uj by Uj, add the 
two equations, and average. This gives, using continuity (eq. (4-21)), 
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Similarly, multiplying equation (4-22) by r, equation (4-23) by Uj, adding the two equations, and 
averaging, give 
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where continuity is again used and buoyancy is neglected. Equations (4-140) and (4-141) are known as 
moment or correlation equations. Setting j = i in equation (4-140) and using continuity, we get, for the 
rate of change of the kinetic energy per unit mass, 
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(4-142) 


ax, ax, ax, ax, 


As in equations (4-25) to (4-27), barred quantities in equations (4-140) to (4-142) which contain lower- 
case letters are turbulent quantities. 

The one-point correlation equations (4-140) and (4-141) give expressions for the rate of change of 
UjUj and Ujr which might be used in conjunction with equations (4-25) to (4-27). But the situation 
with respect to closure is now worse than before. Whereas without equations (4-140) and (4-141) we had 
only to determine u,uj and Ujr , with it we have to determine quantities like UjUjU k , ruju k , <ruj, or, 
a duj/ dxj , (duj/ ax,)(auj/ax,) , o dr/ dxj , and r d 2 Uj/ ax k dx k . One might use equation (4-24) to obtain 
the pressure correlations, but that would only introduce more unknowns. However, equations (4-140) 
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to (4-142) are very useful for studying the processes in turbulence, in that most of the terms have clear 
physical meanings. Moreover, one may be able to calculate terms in those equations from numerical 
solutions of the unaveraged equations. 

An additional one-point equation can be obtained by multiplying equation (4-23) by r, averaging, 
and again neglecting buoyancy. This gives 
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(4-143) 


o t 

where r is the variance of the temperature fluctuation. 

4.3.3. 1 Physical interpretation of terms in one-point equations . — Consider first equations (4-140) 
and (4-142), since physically meaningful interpretations of all the terms in those equations can be given. 

As in the case of equation (3-19) it is helpful, for purposes of interpretation, to multiply the terms in 
equations (4-140) and (4-142) through by p and by a volume element dx 1 dx 2 dx 3 . Then the term on the 

left side of equation (4-140) or (4-142) gives the time rate of change of /?UjUj , or of the kinetic energy 
(OUjUj/2, within the element. This rate of change is contributed to by the terms on the right sides of the 
equations. The first of those terms is equal to the net work done on the element by turbulent stresses 
acting in conjunction with mean-velocity gradients. It is therefore called a turbulence production term; it 
equals the rate of production of /nijUj or of /9 UjUj/ 2 within the volume element by work done on the 
element. A somewhat abbreviated interpretation suggested by the form of the term, which is often given, 
is that it represents work done on the turbulent stress />UjUj by the mean- velocity gradient. 

The next term in each of the equations describes the convection or net flow of turbulence or turbu- 
lent energy into a volume element by the mean velocity U k . It moves the turbulence bodily, rather than 
doing work on it by deforming it, as in the case of the production term. It vanishes when either U k is 
zero (no mean flow) or when the turbulence is homogeneous [ujUj * UjUj(x k )] . In the latter case there is 
no accumulation of turbulence within a volume element, even with a mean flow. 

The next three terms in equation (4-140) and in equation (4-142) also vanish for homogeneous tur- 
bulence. Since they do not contain the mean velocity they do not convect the turbulence bodily or do 
work on it. Therefore, we interpret them as diffusion terms which diffuse net turbulence from one part of 
the turbulent field to another by virtue of its inhomogeneity. The pressure-velocity-gradient terms in 
equation (4-140) drop out of the contracted equation (4-142) because of continuity (eq. (4-21)). There- 
fore, they give zero contribution to the rate of change of the total energy UjUj/ 2 , but they can distribute the 

2 

energy among the three directional components u^/2 (no sum on i). The last term in equations (4-140) 

and (4-142) is the viscous dissipation term, which dissipates turbulence by the presence of fluctuating 
velocity gradients. 

Consider next, equation (4-141) which, on multiplication by pc and the volume element dx x dx 2 dx 3 , 
gives the time rate of change of /?cujr within the element. The first term on the right side gives the 
production of turbulent heat transfer pcujr by the interaction of ru k with the mean velocity gradient and 
by the interaction of UjU k with the mean temperature gradient. The next term gives the net flow of pcujr 
into the volume element by convection. It changes pcujr within the element by moving that quantity 
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bodily. The term vanishes for homogeneous turbulence. The next two terms in equation (4-141) also 
vanish for homogeneous turbulence, but since they do not contain the mean velocity, they do not convect 

P cu i T bodily or produce pcujr . Therefore, as in the analogous terms in equation (4-140), we interpret them 
as diffusion terms which diffuse a turbulence quantity from one part of the turbulence field to another by 
virtue of its inhomogeneity. The next term changes pcx^r by the interaction of the pressure fluctuation 
a with the temperature-fluctuation gradient. The last two terms change pcu^r by molecular action. 

For a = v (for a Prandtl number of 1) they can be interpreted by writing them as v d^Ujf/ dxj dxj 
— 2^ duj/ dx i dr/ dx^. Comparing these terms with analogous terms in equation (4-140), we interpret the 
first as a diffusion term and the second as a dissipation term. 


Finally, consider the equation for the evolution of the temperature variance r 2 (eq. (4-143)). The 
terms on the right side of that equation are similar to terms already considered in the other evolution 
equations. They are, respectively, a production term, a convection term, two diffusion terms which dif- 

fuse r by virtue of the inhomogeneity of the turbulence field, and a dissipation term which dissipates 
2 

r by fluctuating temperature gradients. 


4. 3. 3. 2 Some direct numerical simulations of terms in one-point moment equations . — In section 4.3.3 
we obtained averaged one- point moment (or correlation) equations which give the time-evolution of quantities 
such as UjUj, U|7", and UjUj. Also, in the last section we gave interpretations of the terms in those equa- 
tions. However, because of the closure problem which arises in connection with averaged equations, par- 
ticularly with one-point moment equations, we have not been able to solve those equations or calculate 
terms. Those terms have, however, been calculated by numerical solution of the unaveraged equations 
(eqs. (4-11) to (4-13) (or by eqs. (4-22) to (4-24) together with eqs. (4-25) to (4-27).) 

The cases to be considered here will be concerned with the evolution of the mean turbulent kinetic 
energy according to equation (4-142). The averaged values in that equation will vary only in the direc- 
tion X 2 , one of the directions normal to the mean velocity U^. In that case equation (4-142) becomes 
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Mansour, Kim, and Moin (ref. 49) have calculated the terms in equation (4-144) from a numerical 
solution of the unaveraged Navier-Stokes equations for a steady fully developed turbulent channel flow. 
The boundary condition was u- = 0 at the walls. The results, which correspond to the numerical solu- 
tion plotted in figure 4-2, are shown in figure 4-14. The corresponding turbulence kinetic energy profile is 
plotted in figure 4-15. Since the turbulence is steady-state, the sum of the terms on the right side of 
equation (4-144) is zero. 


The production term — UjU 2 dU 1 /dx 2 , whose form shows that turbulent energy is produced by work 
done on the Reynolds shear stress by the mean velocity gradient, is largest in the region close to the wall and 
peaks at an x 2 of about 12. According to figure 4-2, the mean velocity profile at that point is just 
beginning to deviate appreciably from a laminar (U^ = x 2 ) profile. Note that the production term 
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% u 2 dUi/dx 2 figure 4-14 is positive as it should be, since according to the argument in section 4. 3. 1.1 
and figur e 4-1, when dl^/ dxg is positive UjU 2 is negative, and vise versa. Similarly the dissipation 
term —v du k / dx i du k / dx ( is always negative, because it is the negative of the sum of the squares of 
velocity gradients. The remaining three terms in equation (4-144), and which are plotted in figure 4-14, 
are diffusion terms. They are the kinetic-energy diffusion term du 2 u k u k /2/ dx 2 , the viscous diffusion 
term u d 2 u k u k /2/ dx 2 , and the pressure diffusion term -( davT 2 j dx 2 )/p ■ The plots of those terms show 

that all three of them are positive near x 2 = 0, where a large energy sink occurs because of the boundary 
condition Uj = 0 at x 2 — 0 (see fig. 4-15). Farther away from the wall they become negative. Thus 
they remove turbulent energy from the region where the energy is large and deposit it where the energy is 
smaller. All of the diffusion terms therefore tend to make the turbulence more homogeneous. Far from 
the wall there is a region of small positive du 2 u k u k /2/dx 2 , apparently because the energy in that region 
decreases gradually with increasing x 2 . 

A comparison of the turbulence diffusion processes with spectral transfer processes (to be discussed 
in the next chapter) and the directional-transfer processes arising from the pressure velocity-gradient cor- 
relations (see sect. 4.3.3. 1) is instructive. The spectral-transfer processes remove energy from wave- 
number (or eddy size) regions where the energy is large and deposit it in regions of smaller energy. The 
directional-transfer processes remove energy from large-energy directional components and deposit it in a 
directional component (or components) where the energy is smaller. The turbulence diffusion processes, 
as shown here, remove energy from regions of space where the energy is large and deposit it in regions of 
smaller energy. The spectral transfer, directional transfer, and turbulence diffusion processes tend, 
respectively, to make the turbulence more uniform in wavenumber space, more isotropic, and more homo- 
geneous in physical space. 

Terms in equation (4-144) have also been calculated for a developing turbulent free shear layer, 
where the mean quantities were functions only of x 2 and time (ref. 50). The terms were calculated from 
a numerical solution of the unaveraged equations, where the boundary conditions were periodic. Results 
were similar to those just considered for a fully developed channel flow, insofar as the kinetic-energy and 
pressure-diffusion terms diffused turbulence kinetic energy from regions where the energy was high to 
regions where it was lower. As expected (see the discussion for the channel flow) the production term 
was positive. The viscous dissipation and the viscous diffusion terms were, however, negligibly small, 
unlike those terms for the channel. That was apparently because the presence of walls in the channel 
produced larger mean gradients there than those in the free shear layer. 


4.3.4 Two-Point Correlation Equations 

As shown in section 4.3.3 the use of one-point correlation or moment equations in conjunction with 
the equations for the mean flow causes the number of unknowns to go up faster than the number of avail- 
able equations, so that the closure problem, in effect, gets worse. This leads one to consider the use of 
two-point equations, where one might be able to at least eliminate terms containing instantaneous grad- 
ients, such as o duj/ dxj and du-J dxj duj/ dxj. In that case the use of higher-order equations might be 
expected to increase the number of unknowns at the same rate as the number of available equations. 

Since an assumption for still undetermined higher-order quantities (e.g., triple correlations) might be 
expected to affect the flow to a lesser extent than an assumption for lower-order quantities (e.g., ujUj), 
there might be some advantage in this procedure. Moreover the use of two-point equations enables us to 
consider the important spectral-transfer problem of turbulence, as we will see in the next chapter. 
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Before obtaining the two-point correlation equations, consider the structure of the two-point cor- 
relations. Some of those correlations have been obtained by Kim, Moin, and Moser (ref. 17) in a numer- 
ical solution of the unaveraged Navier-Stokes equations and are plotted in figure 4-16. As the separation 
r k between the defining points increases, the correlation between the velocities at those points goes to 
zero (as r k — > «•). This might be considered a consequence of the mixing property of turbulence; although 
mixing is generally associated with temporal separation (sect. 4.1.1), the effects of temporal and spatial 
separations appear to be similar. In connection with the spanwise separations, note that the correlations 
go negative before becoming asymptotically zero as r 2 increases, in contrast to the correlations with 
point separation in the streamwise direction. More will be said about that phenomenon in the next 
chapter, where it will be predicted theoretically. 


In this section we will consider only two-point equations obtained from the Navier-Stokes equations 
and will neglect heat-transfer and buoyancy effects. Then equation (4-22) for Uj at a point P and a 
similar equation for uj at a point P' can be written as 
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Multiplying equation (4-145) by uj and equation (4-146) by Uj, adding, taking averages, using the fact 
that quantities at one point are independent of the position of the other point, and introducing the new 
variables r k = x k - x k and (x k ) n = nx k + (1 — n)x k (0 < n < 1, see fig. 4-17), results in 
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where the following transformations were used: 
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To obtain the general two-point equations for the pressure-velocity correlations, take the divergence 
of equations (4-145) and (4-146) and use the continuity equation. Thus, from equation (4-146), one obtains 
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Multiplying by Uj, taking averages, and introducing the variables r k and (x k ) , as in equation (4-147), 
gives 
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Similarly, from equation (4-145), 
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The equation for the mean velocity (eq. (4-25)) with buoyancy neglected) is 
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Equations (4-147), (4-149), and (4-150) constitute the two-point correlation equations for inhomo- 
geneous turbulence with mean velocity gradients. A possible solution for these equations is that all the 
correlations are zero. In that case no turbulence exists, and the flow is laminar. We are mainly inter- 
ested here in nontrivial turbulent solutions. 

The terms on the right side of the Poisson equations (4-149) and (4-150) are source terms associ- 
ated with the mean velocity and triple correlations. Most of the terms in equations (4-147) are similar to 
terms in the one-point equation (4-140). Thus equation (4-147) contains turbulence-production, convec- 
tion, viscous, and diffusion terms. The sum of the pressure-velocity correlation terms on the right side of 
the equation transfers energy among directional components. That can be seen by using the relations follow- 
ing equation (4-147) (the first for dn { o' / drj and the second for dcruj/ dr i? and using continuity, since the 
sum of the terms is zero for i = j; the terms give zero contribution to the rate of change of energy UjUj/ 2 , 
but they can transfer energy among directional components. The remaining terms in equation (4-147), 
9(u;Uj'u k ' - UjU k Uj') / dr k and (U k - U k ) du;Uj' / 9r k , are new terms which do not have counterparts in 
the one-point equations. One may interpret them by converting equation (4-147) to spectral form by tak- 
ing its Fourier transform. We will postpone doing that until the next chapter, where the spectral form of 
the equations will be discussed. Note that no averaged quantities containing instantaneous gradients 
appear in those equations, as they do in the one-point equations, so that the number of unknowns, in 
comparison with the number of available equations, is reduced. However, the two-point correlation equa- 
tions, together with the equation for the mean velocity (eq. (4-151)) and the continuity equation (eq. (4-21)) 
still do not form a determinate set for turbulent flow because of the unknown triple correlations. 

If the turbulence is weak (low Reynolds number), and/or if the mean velocity gradients are gen- 
erally very large, it may be possible to neglect terms containing triple correlations. Higher Reynolds 
numbers might be considered by constructing general three-and four-point equations. In each case the set 
of equations would be made determinate by neglecting terms containing the highest-order correlations. 
However, even for the case of a fully developed flow at low Reynolds number, the difficulties of solution 
are extremely great. For that reason modelers of inhomogeneous turbulence have generally preferred to 
use the one-point equations (see e.g., refs. 51 and 52). 4 

Still, it appears that the present scheme should, formally, constitute a solution to the turbulent 
shear-flow problem. But the direct numerical solution of the unaveraged equations, as in references 17, 

49, and 50, may be simpler. One can, however, get an approximate solution for sustained inhomogeneous 
turbulence and study, to some extent, the sustaining mechanism (ref. 53). 

Turbulence is essentially a nonlinear phenomenon.^ The nonlinear character of the two-point correl- 
ation equations (4-147), (4-149), and (4-150), even when the triple correlation terms are neglected, is made 
evident when the mean velocities in those equations are eliminated by introducing equation (4-151) into 
them. A simplified problem is considered, where the mean velocity is in the x 1 -direction and mean quan- 
tities are steady and can change only in the x 2 -direction (ref. 53). Then if the mean velocity is elimi- 
nated by equation (4-151) and the turbulence is weak enough (and/or the mean velocity gradients are 
generally large enough) for triple-correlation terms to be neglected, the set of simultaneous correlation 


4 The two-point equations can, of course be profitably used for statistically homogeneous turbulence, 
as in the next chapter. 

5 The action in turbulence is somewhat similar to that of a clock, a violin bow, or an electronic 
oscillator in that in each of these a steady flow of energy is converted into oscillating energy by a 
nonlinear mechanism. 
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equations which must be considered consists of four nonlinear partial differential equations in four inde- 
pendent and four dependent variables. Although these equations might appear to be too complicated to 
be of use, some results can be obtained by expanding them in power series in each of the space variables 
to obtain algebraic expressions for the correlations, some of which are nonlinear (quadratic). The solu- 
tions so obtained will be accurate only for small values of the space variables. 

The reduction of the partial differential equations to a set of algebraic equations, some of which 
contain quadratic terms, shows the importance of the nonlinear terms in the original equations. The 
quadratic terms in the algebraic equations come, of course, from the nonlinear terms in the differential 
equations. If those terms were absent, the algebraic equations would all be linear and homogeneous (no 
constant terms), since the original equations contained no constant terms. In that case the only solution 
would be that all the turbulence correlations are zero, so that the flow would be laminar. On the other 
hand, in the presence of nonlinear terms the system of algebraic equations is quadratic, so that nonzero 
steady-state turbulence correlations can exist (ref. 53). 

Reference 53 also shows the importance of the pressure-velocity correlations in sustaining turbu- 
lence. The nonlinear production term is absent in the equation for u 2 u 2 ' (i = j = 2 in eq. (4-147)), since 
Ui = in our problem. Thus the energy for sustaining u 2 u 2 must be transferred from the other 

components of UjUj . The equations for the other components do, of course, contain nonlinear production 
terms, and energy can be transferred among the directional components by the pressure-velocity correla- 
tions (see discussion following eq. (4-151)). 

The set of algebraic equations obtained by expanding the correlation equations in power series in 
the independent variables has been reduced to a single quadratic equation 


where 


N 2 = 

C 


N + 


("rt'llt 


N + a(ujU 2 ) 


(4-152) 



N c is the critical value of N (at which UjUj = 0), the subscript 0 designates the point in x 2 , r k space 
about which the series expansions are made, r is the shear stress, which is uniform in the present problem, 
f is a measure of the scale of the turbulence, and a is a constant. Equation (4-152) might be used to 
estimate the variation of u i u 2 with T and {, preferably over a limited range; otherwise a might not 
be constant. For values of N below the critical value the correlation changes sign, according to equa- 
tion (4-152). In that case one should use the no-turbulence solution of equations (4-147) and (4-149) 
to (4-151). Note that the no-turbulence solution is valid when N is above its critical value as well as 
when it is below it. Thus the fluid could be either turbulent or nonturbulent when N is above N c . 

In order to make a rough comparison of our results with experiment we can recast equation (4-152) 
so as to obtain an approximate expression for the friction factor for low Reynolds-number fully developed 
pipe flow. Actually, equation (4-152) is more nearly applicable to Couette flow than to pipe flow; how- 
ever, because of the similarity of trends in pipe and Couette flows and the availability of experimental 
data for transition pipe flow, the comparison will be made with pipe-flow data. The approximate equa- 
tion obtained in reference 53 is 
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f = 16 + BRe 2 Re -1 - BRe~ 3 , 

\ c 


(4-153) 


where f = 2r w /(/>U 2 ) is the friction factor, Re = U a D/z/, U a is the averaged pipe velocity, D is the 

pipe diameter, Re c is the critical Reynolds number, and B is considered constant. For Re = Re c , 
f = 16/Re c , which is the value of f for laminar flow at Re c . 

Equation (4-153) and the solution for laminar flow, together with experimental data from refer- 
ence 54, are plotted in figure 4-18. The transition Reynolds number Re c was taken as 2200, and B 
was set equal to 0.21 xlO 9 . With these values for Re c and B, equation (4-143) is in reasonable agree- 
ment with the data for the turbulent region up to a Reynolds number of about 5500. The deviation of 
the curve from the data at higher Reynolds numbers is probably caused mainly by the neglect of triple 
correlations. The laminar solution is, of course, part of the present solution, since the case of no turbu- 
lence satisfies the correlation equations. For Re < Re c the laminar solution is appropriate because the 
turbulent solution goes below it and thus becomes unphysical. 

4.4 REMARKS 

We began this chapter by considering kinds of averages and their properties, together with some 
notions of randomness from ergodic theory. Next, averaging was applied to the continuum equations 
derived in the last chapter. It was noted that when the equations are averaged, random chaotically fluc- 
tuating variables are replaced by variables which change in a smoother fashion. However, the closure 
problem then arises. Whereas the original unaveraged equations, together with appropriate initial and 
boundary conditions could, in principle, be completely solved, the solution of the averaged equations is 
indeterminate. That is, because of the nonlinearity of the basic equations, there are more unknowns than 
equations in the set of averaged equations. The addition of higher-order moment equations does not alle- 
viate the problem; the number of unknowns goes up as fast as the number of equations when the added 
moment equations are multipoint, whereas if the added equations are single-point, the number of unknowns 
goes up faster. However, the use of those equations enables one to study the processes in turbulence. In 
particular, the numerical solution of the instantaneous (unaveraged) equations enables one to calculate 
terms in the averaged evolution equations. That procedure, in fact, provides a means of closing the sys- 
tem of averaged equations without introducing information which is supplemental to the basic continuum 
equations. 

Because of the engineering importance of obtaining solutions for turbulent flows at high Reynolds 
numbers (where accurate numerical simulations cannot yet be carried out), much effort has been devoted 
to modeling the undetermined terms in the averaged equations. To this end, we obtained some approxi- 
mate solutions for a simple fully developed flow and heat transfer which were, however, based partially on 
exact information. Other approximate solutions and closure schemes, including a frozen stress model for 
highly accelerated flows, were also discussed. The basis for a mixing-length theory and the nature of tur- 
bulent mixing were considered. Finally, the functional forms of the velocity and temperature profiles near 
and away from a wall, and the logarithmic profiles for the velocity and temperature were derived. 

However, the only turbulent solutions considered which could be called completely deductive are 
numerical solutions of the unaveraged equations. As discussed, it is necessary to use supplementary infor- 
mation if the equations are averaged, unless, of course, they are closed by a numerical solution of the ori- 
ginal unaveraged equations. Attempts to obtain analytical solutions of the unaveraged equations have 
not been successful, mostly because of the nonlinearity of those equations. Sensitive dependence of those 
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equations on initial conditions, in fact, appears to preclude an analytical solution. It is mentioned by 
Herring (ref. 55) that the simplest turbulence theory is just the Navier-Stokes equations; since most 
turbulence calculations are numerical anyway, no insight is lost by considering direct integration of the 
Navier-Stokes equations forward in time, starting with some suitable initial data. Because of its rele- 
vance and importance in turbulence research, we give here a few additional remarks on direct numerical 
solution. 

Although numerical solution (or numerical simulation) bears some similarity to experiment there is 
an important difference. The former uses directly, and attempts to solve, a given set of mathematical 
equations (say the Navier-Stokes equations). The latter ordinarily does not, although both methods may 
arrive at the same description of nature. A numerical solution is, in fact, much closer to a theoretical 
solution than to experiment. A numerical solution differs from a theoretical one only in the tools that are 
used; numerical methods and high-speed computers may take place of, or supplement, analytical tools. A 
significant difference is not, as sometimes stated, that the numerical solution is discrete; as many points 
as desired can be calculated if a suitable interpolation formula is provided. 

Because of the small scale of some of the turbulent eddies, accuracy is a problem in the numerical 
study of turbulence. Smaller and smaller eddies are generated as the Reynolds number (strength) of the 
turbulence increases. One can always pick a Reynolds number large enough that the results will be quan- 
titatively inaccurate, no matter how fine the numerical mesh. On the other hand, turbulence at any fixed 
Reynolds number could, in principle, be calculated, given the availability of a powerful enough computer. 
Note that there does not seem to be a qualitative difference (no bifurcations) between low and high 
Reynolds-number turbulence; in the latter the turbulent energy is just spread out over a (much) wider 
range of eddy sizes. The mathematical or computational methods used may of course, be different for 
higher Reynolds numbers. 

Although we started this section with a discussion of averages and averaged equations, it has been 
appropriate to close it with a discussion of numerical solutions of the original unaveraged equations, those 
being the only completely deductive turbulent solutions considered thus far (see, e.g., refs. 17 and 49). 
(Deductive refers here to solutions of the continuum equations which do not require an input of supple- 
mentary information.) Note that, as mentioned earlier, a deductive solution of an unclosed averaged 
equation (e.g., the evolution equation for the turbulent energy) can be obtained by using it in conjunction 
with unaveraged equations, the terms in the averaged equation being calculated from the numerical solu- 
tion of the determinate unaveraged equations. As a final remark, it seems unlikely that a turbulent 
analytical solution of the unaveraged equations could ever be obtained, because the solution would have 
to be sensitively dependent on initial conditions (see chapter VI). 

In order to continue with the analysis of turbulence, it will be necessary, or at least convenient, to 
introduce Fourier analysis. That will be done in the next chapter. 
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Figure 4-1. — Ske tch illustrating development of turbulent shear 
stress — pu,u 2 or heat transfer pc -fu 2 in presence of a mean 
velocity or temperature gradient. Regardless of whether the 
velocity fluctuation u 2 is positive or negative, the product "ilj7u 2 
(or tOj) is negative (see section 4.3.1 .1). Also illustrated (in the 
lower part of the figure) is a mixing-length theory, where e 2 Is 
the mixing length and x 2 0 is at the virtual origin of an eddy (see 
section 4.3.2.2). 
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Figure 4-2. — Velocity distribution for turbulent flow through a channel or pipe. 
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Figure 4-3. — Comparison of experimental fully developed turbulent friction factors for pipe flow with those calculated from eqs. (4-106), 
(4-112) and (4-113). 






.001 




Transfer 

Ref. 


O Ethylene glycol and water 

Heat 

20 

J— 

□ Water 

Heat 

21 

E 

O Water 

Heat 

32 

■ 

> Air 

Heat 

23 


V Sodium hydroxide 

Heat 

24 


D Sodium hydroxide 

Heat 

25 


D Zinc in mercury 

Mass 

26 

— 

O Salts in water 

Mass 

27 

£ 

o Vaporization 

Mass 

28-30 

- 

D Diffusion-controlled electrode Mass 

31 


.0001 t 


(Each type of symbol represents a mean through data) 

Calculated from eqs. (4-1 06), (4-1 1 6), 

and (4-104) (c/e h = 1) 


.00001 


L-l. 1 I lllll L-i I ill 111 I i I i I ilil 

(a) Re = 2 (x 2 ) c = 1 0 000. 


r T9 

I 

i 

Y 

o 

ii 

55 


.001 




Transfer 

Ref. 

O 

Ethylene glycol and water 

Heat 

20 

□ 

Water 

Heat 

21 

o 

Water 

Heat 

32 

> 

Air 

Heat 

23 

7 

Sodium hydroxide 

Heat 

24 

O 

Salts in water 

Mass 

27 

o 

Vaporization 

Mass 

28-30 

D 

Diffusion-controlled electrode 

Mass 

31 


(Each type of symbol represents a mean through data) 


.00001 



.0001 u 


Calculated from eqs. (4-106), (4-116), 
and (4-104) (€/e h = 1) 



(b) Re = 2 (x 2 ) 0 U^v 

= 25 000. 



.01 




Transfer 

Ref. 


“ 

O 

Ethylene glycol and water 

Heat 

20 



□ 

Water 

Heat 

21 


- 

< 

Water 

Heat 

22 



^ o 

Water 

Heat 

32 




Air 

Heat 

23 

.001 



° 

Vaporization 

Mass 

28-30 


~ 


Aniline 

Heat 

22 


7 

° 

n-Butyl alcohol 

Heat 

22 


.0001 


(Each type of symbol represents a mean through data) 

Calculated from eqs. (4-106), (4-116), 

___ and (4-104) (t/t h = 1) 

Calculated using e/e h = 0.7 


.00001 


.1-1 1 I M ill 


i.lilil 


100 


Pr = v/a 


1 000 


10 000 


(c) Re = 2 (x 2 ) c Ua/v = 50 000. 

Figure 4-4.— Comparison of calculated Stanton numbers with experiment. (For the mass-transfer data the tem- 
peratures T w and T a are replaced by concentrations , and the molecular thermal diffusivity a is replaced by a 
diffusivity for mass transfer.) 
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(a) Blackwelder & Kovasznay (ref.40). 



(b) Patel & Head (ref. 44) (Note shifted vertical scales). 

Figure 4-5.— Predicted early and intermediate development of mean velocity profile 
in a turbulent boundary layer with severe favorable pressure gradients and com- 
parison with experiments. Symbols indicate experimental data points. 
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theory equation (4-135) 
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Figure 4-6. -Semi-logarithmic law-of-the-wall plot of 
theoretical velocity profiles for severe favorable pressure 
gradients and comparison with experiment of Patel & 
Head (ref. 44). (Note shifted vertical scales.) Symbols 
indicate experimental data points. 
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Figure 4-7. —Effect of neglecting Reynolds shear stress on 
theoretical law-of-the-wall plot for experimental conditions 
of figure 4-5. (Noted shifted vertical scales.) 
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Figure 4-8. — Contribution of various terms in eq. (4-1 28) to streamwise rate of 
change off dimensionless kinetic energy of mean flow for initial conditions of 
figure 4-6. Terms in eq. (4-128) nondimensionalized by Uoand p. U„/Uo =1.14. 

Note that transverse component of kinetic energy is negligible for a boundary 
layer. 
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Figure 4-9. — Effect of mass injection normal to wall on theoretical velocity profiles for severe favorable pressure 
gradients. (Initial conditions and pressure gradients correspond to figure 4-6) 
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Figure 4-10.— Theoretical velocity profiles for severe adverse pressure gradients and comparison with experiment 
of Kline et al. (ref. 45). (Note shifted vertical scales). 


58 




1 


2 


3 


4 


5X10-5 


0 


U 0 (x, - xJ/m 


Figure 4-1 1 .—Comparison of streamwise evolution of theoretical Stanton number St and 
skin-friction coefficient Cfin moderately short highly accelerated turbulent boundary layers. 
Xo is value of x, at initial station. U, is value of at initial station. 





Figure 4-12. — Comparison of evolution of theoretical velocity and temperature distributions in 
moderately short highly accelerated turbulent boundary layers. 
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Figure 4-1 3. — Comparison of theory and experiment for streamwise evolution of Stanton number in moderately 
short highly accelerated turbulent boundary layers. Xq is value of x, at initial station, fv, is maximum value of 
(v/pUoo 3 ) dP/dx, for a particular run. 
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Figure 4-14.— Plot of terms in one-point equation for evolution of kinetic energy (eq. (4-144)) for fully 
developed channel flow. Results obtained by numerical solution of unaveraged Navier-Stokes eaua- 
tions (Mansour, Kim, and Moin ref. 49). Terms nondimensionalized by/f^Tand v. Channel Reynolds 

number (based on centerline U, and channel half-width) is 3200. 
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Figure 4-15.— Plot of turbulence kinetic energy for fully developed channel flow in figures 4-14 and 4-2. Results obtained by numerical 
solution of unaveraged Navier-Stokes equations (Kim, Moin, and Moser, ref. 17). 
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Figure 4-16 — Variation of two-point correlation coefficients with point seperation as obtained by Kim, Moin, and Moser (ref. 17) from a 
direct numerical simulation of unaveraged Navier-Stokes equations. x^(x^ c = 0.030, x 2 + = 5.39. 
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Figure 4-17. — Vector configuration for two-point correlation 
equations. 



Figure 4-1 8.— Comparison of eq. (4-1 53) and solution for laminar 
flow, with pipe-flow data in turbulent transition and laminar 
regions. 
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